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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 


A  multiple  foil  activation  iterative  method  has  been  developed  to  determine 
neutron  flux  spectra.  A  computer  code  (SAND  II)  has  been  written  which  pro¬ 
vides  a  "best  fit"  neutron  differential  flux  spectrum  for  a  given  input  set 
of  infinitely  dilute  foil  activities.  The  results  of  experimental  and  ana¬ 
lytical  studies  for  a  wide  variety  of  neutron  environments  strongly  suggest 
that,  in  addition  to  determining  neutron  flux  spectra,  the  multiple  foil 
activation  iterative  method  can  be  helpful  in  (1)  the  validation  and  improve¬ 
ment  of  calculntional  techniques  used  to  predict  neutron  flux  spectra,  (2) 
upgrading  of  confidence  in  neutron  spectrometer  measurements,  (3)  determining 
the  reliability  of  foil  activity  measurement  methods,  (4)  assessing  material 
scattering/absorption  effects,  and  (5)  examining  current  foil  detector  cross 
section  evaluations  to  provide  guidance  for  re-evaluation  for  these  data,  and 
eventual  "unfolding"  of  the  absolute  differential  form  of  cross  sections  for 
any  foil  reactions  producing  a  detectable  product . 
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INTRODUCTION 

The  authors  have  previously  reported  an  iterative  perturbation  method  for 
determining  neutron  flux  spectra  by  multiple  foil  activation  (Refs.  1,2,  3).  This 

method  required  certain  subjective  decisions  on  the  part  of  the  user  to  gener- 

>;< 

ate  each  perturbation.  A  computer  code  SAND  II  has  since  been  developed  to 
determine  spectra  by  a  fully  automated  iterative  method,  described  in  Sec¬ 
tion  III.  Greer  and  Walker  (Ref.  4)  have  studied  this  same  problem  from  a 
somewhat  different  point  of  view  and  have  developed  a  code  (SPECTRA)  which 
gives  integral  results  comparable  to  those  of  SAND  II.  Both  of  these  methods 
have  been  subjected  to  extensive  analytical  and  experimental  testing,  the  re¬ 
sults  of  which  (for  the  SAND  II  code)  are  discussed  in  Sections  IV  and  V.  Eoth 
the  immediate  and  long-range  implications  of  this  work  for  problems  related  to 
neutron  spectral  definition  are  considered  in  Section  II. 

The  SAND  II  code  provides  a  "beat  fit"  neutron  flux  spectrum  for  a  given 
input  set  of  infinitely  dilute^  foil^  activities.  The  calculational  procedure  con¬ 
sists  of  the  selection  of  a  known  flux  spectrum  form  to  serve  as  the  initial 
approximation  to  the  solution,  and  subsequent  iteration  to  a  form  acceptable  as 
an  appropriate  solution,  The  solution  is  specified  either  as  time-integrated 
flux  (fluence)  for  a  pulsed  environment  or  as  flux  for  a  steady-state  type  neu¬ 
tron  environment. 

The  results  of  the  analytical  and  experimental  studies  presented  in  Sec¬ 
tions  IV  and  V  indicate  that  the  multiple  foil  iterative  method  of  neutron  l'lux 
spectral  determination  can  be  expected  to  give  integral  neutron  flux  results 
over  the  energy  range  from  1Q~^  Mev  to  18  Mev  that  are  accurate  to  within 
±10%  to  30%  at  all  energies,  if  the  errors  in  the  foil  reaction  cross  sections 
and  measured  activities  are  within  similar  .limits.  (The  problem  of  errors  is 


^Spectrum  Analysis  by  Neutron  Detectors  II. 

tlnfinitely  dilute  means  that  the  density  of  target  nuclei  is  small  enough  to 
preclude  self-absorption  and  scattering  effects. 

§The  word  "foils"  as  used  in  this  report  is  meant  to  include  all  physical  forms 
in  which  target  nuclei  of  any  particular  element  are  irradiated,  foils,  pellets, 
wires,  powder,  etc. 
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given  consideration  throughout  this  Volume,  and  particularly  in  Section  VI.) 
Differential  flux  spectra  are  also  obtained  which  reflect  true  structure  as  well 
as  fluctuations  resulting  from  experimental  errors.  The  method  is  such  that 
the  use  of  a  large  number  of  foils  with  overlapping  energy  regions  of  sensitivity, 
and  subsequent  examination  of  the  solution  differential  spectral  structure,  can 
be  used  to  distinguish  the  true  structure  from  that  which  may  be  caused  by  cross 
section  errors,  activity  errors,  and  absorption  effects. 

SECTION  II 

FLUX  SPECTRAL  DEFINITION 
A,  Review  and  Objectives 

Radioactivation  techniques  are  in  common  use  for  measuring  neutron  flux, 

If  the  neutron  spectrum  can  be  well  represented  as  a  fission  spectrum,  straight¬ 
forward  procedures  are  available  for  the  interpretation  of  the  experimental  data 
(Ref,  5).  If  the  form  of  the  neutron  spectrum  is  not  known,  or  known  only  ap¬ 
proximately,  the  problem  of  interpretation  is  more  difficult  and  at  present  is  the 
subject  of  considerable  review  and  investigation  (Refs,  6,7), 

For  fast  neutron  flux  spectra  (i,  e,  ,  E  >  0,1  Mev),  methods  that  involve  the 
use  of  a  set  of  materials  with  different  neutron  energy  sensitivity  limits  have 
been  proposed  and  studied  quite  extensively  (Refs.  6,8,9,10,  11,12),  Difficulties 
in  interpreting  the  experimental  data  arise,  however,  caused  by  the  propagation 
of  experimental  errors  in  the  system  of  activation  integral  equations.  Those 
errors,  associated  with  uncertainties  in  available  cross  section  data  and  in  the 
measured  values  of  saturated  activities,  place  limitations  on  the  mathematical 
techniques  that  can  successfully  be  used  to  determine  neutron  spectra  from  the 
foil  activation  '  tu,  (In  the  present  work,  realistic  UrnitH  have  been  assigned 
to  these  e"r<  5th  for  the  experimental  and  analytical  studios,) 

The  cur  re  date  of  technology  of  neutron  flux  spectral  definition  includes 
problems  in  f  .  1  areas,  affecting  both  calculationul  and  experimental  tech¬ 

niques.  The  ireas  are 


1,  Reaction  Cross  Sections 


Current  cross  section  evaluations  for  neutron  interactions  used  in  spec  ¬ 
tral  calculations  and  measurements  by  both  spectrometers  and  the  foil  activation 
technique  contain  uncertainties  and  inaccuracies  in  both  absolute  magnitude-  and 
differential  form. 

2.  Activity  Measurements 

Measurements  of  foil  activities  are  hampered  by  incomplete  knowledge 
of  decay  schemes,  interfering  impurities,  uncertainties  in  half-life,  nonavail¬ 
ability  of  standards  for  short  half-life  products,  etc, 

3.  Absorption  Effects 

Absorption  effects  may  be  manifest  as  apparent  errors  in  activities  for 
the  foil  technique  and  as  a  decrease  in  energy  resolution  for  spectrometers, 

4.  Spectrometry 

Accuracy  of  spectrometer  measurements  iB  affected  by  such  factors  as 

3 

uncertainties  in  cross  section  (such  as  He  ),  sensitivity  of  electronics ,  calibra¬ 
tion  of  efficiency,  and  sensitivity  to  secondary  radiation,  such  as  gamma  rays 
in  a  neutron  environment. 

5.  Calculation 

There  is  a  need  to  prove  the  validity  of  calculational  techniques  for  neu¬ 
tron  environments,  in  general,  over  the  whole  neutron  energy  range.  The  iter¬ 
ative  foil  activation  method  Is  uniquely  suited  to  the  problem  of  the  validation  of 
analytical  methods  that  have  been  developed  for  the  characterization  of  neutron 
flux  spectra,  by  allowing  a  direct  comparison  with  experiment  (see  Section  V). 
Such  validation  is  necessary  to  give  confidence  in  the  analytical  calculations  that 
will  be  required  for  locations  in  which  the  placement  of  foils  or  spectrometers 
is  not  possible  or  practicable,  due  to  such  considerations  as  high  temperature, 
accessibility,  or  space  limitation. 

The  authors  strongly  believe  that  through  a  marriage  of  the  results  of 
foil  activation  measurements,  spectrometry  and  calculational  techniques,  signifi¬ 
cant  progress  can  be  made  in  resolving  some  of  the  current  problems  in  tin- 
definition  of  neutron  environments  and  subsequent  studies  of  neutron  material 
damage  effects. 
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The  immediate  objectives  of  the  present  work  relating  to  neutron  spec¬ 
tral  definition  have  been  ( 1 )  to  develop  a  computer  code  (SAND  II)  for  determina¬ 
tion  of  neutron  flux  spectra  from  foil  activity  measurements  by  iterative  solution 
of  a  set  of  activity  integral  equations,  (2)  to  increase  knowledge  of  neutron  flux 
spectral  forms  for  a  wide  range  of  neutron  environments,  (3)  to  examine  current 
cross  section  data  evaluations  for  foil  reactions  used  in  neutron  spectral  deter¬ 
mination  and  to  provide  guidance  for  re-evaluation  for  these  data,  (4)  to  deter¬ 
mine  the  reliability  of  current  foil  activity  measurement  methods,  and  (5)  to 
establish  direction  for  future  application  of  the  foil  activation  technique, 

Long-range  objectives,  toward  which  this  program  has  been  the  first 
step,  include  the  improvement  of  neutron  flux  spectral  dosimetry  by  foil  activa¬ 
tion,  upgrading  of  confidence  in  neutron  spectrometer  measurements,  validation 
and  improvement  of  calculational  techniques  for  flux  spectra,  and  "unfolding"  of 
the  absolute  differential  form  of  cross  sections  for  any  foil  reactions  producing 
a  detectable  product, 

B,  Foil  Activation  Equations 

The  foil  method  of  neutron  spectral  measurement  depends  upon  neutron 
activation  of  a  set  of  foils  and  the  subsequent  determination  of  the  radioactive 
decay  rates  (activities)  of  the  foils.  A  general  mathematical  representation  of 
this  process  must  consider  the  differential  flux  as  a  function  of  both  energy  and 
time  [pOM)], 

Let  "n"  be  the  number  of  neutron  reaction  product  atoms  present  in  a  given 
target  foil  material  at  any  time  t,  and  "m"  be  the  number  of  target  atoms,  The 
net  time-rate  of  production  of  reaction  product  atoms  Is 

=  mEFmn(t)  *<*)  ‘  n0rn(O*(t)  -  n*n  (1) 

where  An  is  the  radioactive  decay  constant  of  the  neutron  reaction  product, 

~d  (t)  is  the  spectrum-averaged  cross  section  for  the  reaction  producing  the  " 
inn  in 

atoms,  8^(t)  is  the  total  spectrum-averaged  cross  section  for  all  neutron  reac¬ 
tions  that  remove  any  of  the  "n"  atoms  and  <J>(t)  (often  called  "integral  flux'  )  is 
the  time-differential  flux  defined  by  the  equation 
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n 


(2) 


r  2 

4>{l)=i  <p(E  t)  dE,  (n/'cm-sec) 

■'o 

The  time-rate  of  change  of  target  atoms  (burn-up)  is 

T3T  -  till  (3) 

where  0^(t)  is  the  total  spectrum-averaged  cross  section  for  all  neutron  reac¬ 
tions  that  remove  any  of  the  "m"  atoms.  Equation  (3)  has  the  solution 

-/'‘a L(t')  4(t')  dt' 

J0  ,  A  V 

m  a  mQ  e  (4) 

The  integration  constant  mQ  is  the  initial  number  of  target  atoms  or  the  number 
of  target  atoms  at  time  t  =  0. 

In  Eqs.  (1),  (3),  and  (4),  the  spectrum-averaged  cross  sections  are  defined 
by  equations  of  the  form 


where  <JX(E)  is  the  energy-dependent  cross  section. 

Combination  of  Eqs.  (1)  and  (4)  gives 

&  =  "oW*1  *(t»  0  -”(V  VH  161 

Equation  (6)  is  an  ordinary  first  order  first  degree  linear  differential  equation 
which  has  the  solution 
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(7) 


n  =  in  e 
o 


|‘V  -f'a n(t')dt'i  ft 

o 

X  t'  +  /"V  (t")  -ff  (t")|  dt" 
n  Jin  m  ] 


dt' 


Under  the  restriction  that  there  is  no  significant  burn-up  of  the  target  or  product 
atoms  (that  is,  both  CT^ft )  4>(t)  and  0^(t)  $(t)  «  X^),  Eq.  (7)  reduces  to 


-X  t  -t  X  t' 

n  =  mo  e  n  J  *<t')  e  n  dt' 


(8) 


Two  different  extremes  in  the  flux  time -behavior  require  consideration: 
(1)  steady- state,  and  (2)  pulsed.  In  the  steady-state  case, 


$<t>  =  $ 


(9) 


and 


SLJt) 

mn 


mn 


(10) 


That  is,  both  flux  and  spectrum-averaged  cross  section  are  constant  with  time 
[see  Eq.  (5)J,  Therefore,  Eq.  (8)  may  be  rewritten  as  follows: 


-X  t  rt  x  t' 

n  ~  m  <7  $  e  /  e  dt' 

o  mn  .L 


(11) 


which,  upon  integration,  gives 


n1  iwn  n  ^ 

n  =  “x —  * 

n 


(id) 


The  radioactive  decay  rate,  or  activity,  of  the  product  material  at  time  t  is 
given  by 

a(t)  =  Ann  =  mo?Tmn  *(l  -  e  n  )  (13) 

The  saturated  activity  A  (the  activity  for  an  infinitely  long  irradiation)  may  be 
expressed  by  rewriting  Eq.  (13)  as  follows: 


A  =  a^U-rr  =  mV  « 
o  m  n 


RJ 

Equations  (5),  (9),  (10),  and  (14)  may  be  combined  to  give 

A  -  mo|°°(Tmn(E)<p(E)  dE 


(14) 


(15) 


where 


<p(E)  =  <p( E,  t) 


(16) 


since  flux  is  constant  with  time. 

If  a  set  of  Eq.  (15)  is  obtained  by  irradiation  of  several  foils  with  different 
reaction  energy  responses,  a  matrix  is  provided  from  which  <p(E)  may  be  deter¬ 
mined.  The  solution  of  this  matrix  is  considered  in  Section  III. 

In  the  case  of  a  pulsed  environment,  the  time  behavior  of  4>(t)  is  in  general 
not  constant,  and  may  not  be  measurable.  Equations  (9)  and  (10)  do  not  apply, 
therefore,  and  Eq.  (8)  cannot  be  simplified  to  Eq,  (11).  For  a  pulse,  however, 
Ant  «  1.  Therefore,  both  exponential  terms  in  Eq.  (8)  may  be  approximated  by 
unity,  and  Eq,  (8)  reduces  to 


(17) 
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Activity  is  given  by 


a(t)  =  X  n=X  m  f  a  ( t ' )  $(t‘’)  dt'  (18) 

n  n  ol  mn  ' 

where,  as  before,  a. ( t)  is  the  activity  measured  at  time  t  (in  this  case  at  the  end 
of  the  neutron  burst  or  pulse),  Equations  (5)  and  (18)  may  be  combined  and  re¬ 
written  as  follows: 


a(t)  =  /V  (£)/%{£,  t')dt'  dE  (19) 

"o  Jo 

Comparison  of  Eq,  (19)  for  the  pulsed  case  with  Eq.  (15)  for  the  steady- 
state  case  shows  that  the  same  numerical  procedure  can  be  used  to  treat  both 
cases.  In  the  steady-state  case,  saturated  activities  are  calculated  from  a  dif¬ 
ferential  flux;  in  the  pulsed  case,  end-of-burst  activities  are  calculated  from  a 
differential  fluence.  (For  the  pulsed  case  there  is  no  requirement  that  <p(E.  t) 
be  separable.) 

The  SAND  II  code  has  the  option  of  obtaining  a  spectral  solution  as  differ¬ 
ential  flux  for  the  steady-state  case,  or  as  differential  fluence  for  the  pulsed 
case  (see  Vol.  II,  Section  II-A-2),  Except  for  the  multiplication  of  the  calculated 
integral  by  the  decay  constant  for  the  pulsed  case  [see  Eq.  ( 1 9 ) J ,  and  the  ex¬ 
trapolation  to  saturated  input  activities  for  the  steady-state  case,  the  problems 
of  obtaining jf1  <p(E,t*)  dt*  andtp(E)  respectively,  are  the  same. 

The  discussion  of  the  SAND  II  iterative  method  in  Section  III  is  based  on  the 
equations  for  the  steady-state  case, 

SECTION  III 

SAND  II  ITERATIVE  METHOD 


A.  General  Description 

The  SAND  II  code  is  designed  to  provide  a  "best  fit"  neutron  flux  spectrum 
for  a  given  input  set  of  infinitely  dilute  foil  activities.  As  explained  in  Section 
II-B,  the  input  activities  are  extrapolated  to  saturation  in  the  case  of  a  steady- 
state  environment,  i.  e,  ,  reactors,  and  are  actual  time -of- removal  activities  in 


-8- 


the  case  of  a  pulsed  environment;  the  solution  spectrum  is  differential  flux  in 
the  steady-state  case  and  differential  fluence  in  the  pulsed  case. 

The  code  uses  a  discrete  interval  description  of  the  neutron  energy  range, 
rather  than  a  continuous  representation  by  analytical  functions;  the  energy  range 
between  10'^  and  18  Mev  is  represented  in  620  intervals  (45  per  decade  up  to 
1  Mev,  and  170  between  1  and  18  Mev).  The  solution  spectrum  is  thus  presented 
in  tabular  form  at  621  points.  The  problem  is  essentially  to  solve  for  621  un¬ 
knowns  (solution  differential  flux  values)  in  a  system  of  n  linear  activity  equa¬ 
tions,  where  n  is  the  number  of  foils  used. 

The  calculational  procedure  consists  of  selection  of  an  initial  approximation 
spectral  form  and  iteration  from  that  approximation  to  a  form  acceptable  as  an 
appropriate  solution.  The  solution  is,  of  course,  not  unique,  since  the  number 
of  equations  is  much  smaller  than  the  number  of  unknowns  (n  «  621);  appropri¬ 
ateness  of  the  solution  therefore  depends  on  a  suitable  choice  of  the  initial  ap¬ 
proximation,  based  on  all  available  physical  information  in  any  given  case.  A 
mathematical  statement  of  the  problem  and  method  of  solution  is  presented  in 
Section  III-B.  Briefly,  the  iterative  procedure  consists  of  the  following  steps: 

1.  Activities  are  calculated,  for  each  foil  used,  based  on  the  current 
iterative  spectrum  and  the  evaluated  reaction  cross  section  library 
which  is  part  of  the  code.  (See  Volume  II  for  additional  information, 
and  Volume  III  for  graphical  and  tabular  presentations  of  the  evalu¬ 
ated  cross  sections  as  a  function  of  energy  for  the  29  foil  reactions 
in  current  use. ) 

2.  These  calculated  activities  are  compared  with  the  measured  activities 
to  provide  a  correction  factor  for  each  foil. 

3.  A  weighting  function  (of  energy)  is  obtained,  for  each  foil,  based  on  the 
sensitivity  function  (differential  cross  section  multiplied  by  differential 
flux)  for  that  foil  calculated  at  the  current  iteration. 

4.  The  foil  weighting  functions  are  applied  to  an  averaging  procedure, 
to  obtain  an  average  correction  factor  at  each  energy,  based  on  the 
comparison  of  measured  to  calculated  activity  for  each  foil  and  on  the 
relative  contribution  of  the  flux  at  the  given  energy  to  the  activity  of 
that  foil. 


> 
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5.  The  average  correction  factors  are  then  applied  to  the  current  iterative 
flux  value  at  each  energy,  to  obtain  the  next  iterative  flux  spectrum. 

The  criterion  currently  incorporated  into  the  code  for  achievement  of  an 
acceptable  solution  is  based  on  a  comparison  of  successive  iterative  differential 
flux  values;  a  solution  is  considered  to  have  been  achieved  when  the  difference 
between  two  successive  values  is  smaller  than  a  specified  (as  input)  percent  at 
all  621  energy  points.  (See  Section  IV-B-2.  ) 

SAND  II  includes  the  option  of  selecting  an  initial  approximation  spectral 
form  from  among  a  library  of  reference  forms,  by  searching  for  the  spectrum 
for  which  the  set  of  calculated  activities  for  the  foils  used  is  most  nearly  con¬ 
sistent  with  the  input  (measured)  activities.  Since  no  physical  information  is  in¬ 
corporated  into  this  selection  method,  it  is  emphasized  that  the  reference  spec¬ 
trum  library  is  meant  to  be  available  for  use  only  as  a  last  resort,  in  cases  fur 
which  so  little  information  is  available  that  no  meaningful  selection  of  an  initial 
approximation  can  be  made  on  the  basis  of  physical  knowledge  of  the  measured 
environment.  The  reference  spectrum  library  is  also  an  important  convenience 
in  certain  analytical  studies  aimed  at  understanding  the  solution  behavior  of  the 
iterative  method,  in  calculating  spectrum-averaged  cross  sections,  and  in  cal¬ 
culating  anticipated  activities  and  foil  energy  sensitivity  limits  for  planned  ir¬ 
radiations.  (See  Volume  II  for  additional  information,  and  Volume  IV  for  graphi¬ 
cal  and  tabular  presentations  of  differential  flux  for  the  59  reference  spectra  in 
current  use,) 

B,  Mathematical  Representation 

The  mathematical  procedure  for  the  present  method  involves  iterative 
spectral  perturbation  to  obtain  a  "best  fit"  solution  for  a  system  of  simultaneous 
activation  integral  equations  obtained  from  Eq.  (15)  with  certain  notational  addi¬ 
tions  and  modifications,  defined  as  follows; 

t  h 

A.  s  measured  activity  for  the  i  foil  reaction  (extrapolated  to 

saturation  and  infinite  dilution); 

A[kj  £  calculated  activity  for  the  i1'"1  foil  reaction,  based  on  the  k**1 

iterative  spectrum; 

C>'R^(E)  £  ktR  iterative  differential  flux; 
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tr 


E. 

J 

^  rvi 


A.L 


[k] 


i.j 


<yrm 


N 

i 

Xt 

j 

i 

r 

k 


j  +  1 


m  energy  of  the  energy  point; 

1  ^ 

•  integral  flux  in  the  j  *'  energy  interval,  between  E.  and  E 
t  H 

25  i  foil  reaction  dross  section  for  a  specified  (n,y),  (n,p).  (n,a). 
(n.n^,  (n,f),  or  (n,  2n)  interaction; 

■  that  portion  of  aP^  contributed  by  neutrons  in  the  j'*1  energy 
interval,  between  and  E^jj 

tK 

■  r  foil  cover  removal  cross  section  for  cadmium,  boron  10 
or  gold  covered  foils; 

til 

■  r  foil  cover  nuclei  density; 
til 

■  r  foil  cover  thickness; 

B  1,  2, . m  (energy  interval  index;  currently,  m  =  620) 

■1,2 . n  (foil  index); 

■  1,  2,  or  3  (foil  cover  index); 

*  1,  2,  3,  ,  (iteration  index); 

(Currently,  Ej  =  10"*®  Mev  and  +  j  =  18  Mev.) 

Each  subscripted  equation  below  is  meant  to  represent  a  set  of  equations  on 
the  subscript(s)  as  defined  above. 

fk] 

For  bare,  cadmium-,  boron  10-,  and  gold-covered  foil  detectors,  A,  , 

M  J 

may  be  given  (when  Ej  j  -  Ej  is  small)  by 


1 


1)J  L 


Ej+i 


ej 


rkl  3  -NrX 
O.(E)  <pLkJ(E)  n  e  r  ,J 
1  r  =  1 


dE 


(20) 


l  If 


where  V  .  is  an  appropriate  average  value  for  cr  (E)  over  the  j1*1  energy  interval, 
r ,  J  r 

An  adequate  approximation  for  &  .  for  the  energy  intervals  in  the  SAND  II  code 

r,  j 

has  been  found  to  be 
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(T  * 

r.J 


/*Ej+^ 

/  cr 

Je.  1 
_ 1 _ 

r 


(E)  riE 


This  approximation  assumes  that  the  foil  cover  interval-averaged  cross  section 
is  independent  of  "k"  and  <p^(E):  errors  introduced  by  this  assumption  src  neg¬ 
ligible  compared  to  those  introduced  by  neglecting  the  complicated  dependence 
of  the  scattering  effect  on  the  geometry  of  each  particular  case,  which  generally 
requires  error-introducing  simplifications.  These  errors  are  in  turn  small,  in 
most  practical  cases,  because  of  the  relative  magnitudes  of  the  attentuation 
effects. 

For  purposes  of  calculation,  it  is  convenient  to  define  the  foil  reaction 
interval-averaged  cross  section: 


/  jH  ikj 

/  a(E)  <PLKj(E)  dE 

Je j  _ 

/  dE 

Je. 


Combining  Eqs.  (20)  and  (22)  gives 


EJ  L' 


p,,kJ<E)  n 

r  =  1 


-n  x  tr 

e  r  rr’JdE 


Because  of  the  fine  structure  of  the  evaluated  cross  section  CJ.(E)  for  many 

jk:  1 

of  the  reactions  used,  the  calculation  of  OT  r  by  Eqs,  (22)  requires,  for  many  i, 
a  much  finer  energy  interval  model  than  is  used  by  the  SAND  II  code,  If,  for 
caleulational  facility,  the  approximation  is  made  that 

yJki(E)  =  for  E.  <  E  <  E.  +  1  (24) 


then,  from  Eqs,  (22)  the  foil  interval-averaged  cross  section  0.  .  may  he  eaten- 

J 

lated  by 


L 


Ej+1 


o\(E)  dE 


U. 


i.J 


fE  j  +  1 

/  dE 


:?;) 


This  calculation  is  k-independent  and  does  not,  therefore,  require  recalcu¬ 
lation  at  each  iteration.  In  fact,  £r.  .  (as  well  as  .)  is  calculated  only  once,  in 
j  separate  code  (CSTAPE)  which  performs  exact  linear  interpolation  integration 
over  ff.(E)  and  <y  (E)  input  data  points  for  all  i  (and  r)  and  j.  If  integral  flux  in  the 

th  1  T 

j  energy  interval  is  defined  by 


rE  j+1 


(E)  dE 


(26) 


equations  (23),  (25),  and  (26)  may  be  combined  to  give 


A^kj  is  now  given  by 


A“  ■  S  j  j,  ;NrXr°’r,J 


,  m 

AW  -  V  AW 

1  t—d  1 J 

j=l 


(27) 


(28 ) 


For  describing  the  SAND  II  iterative  procedure,  the  activity-weighting 

Fki 

function  W.  .  is  introduced  and  defined  using  Eqs.  (27)  and  (28): 


*•  A.^  ,VA.Ck]  ,  j  =  2, 

i,J  2^  l, j  i,J-l /  l  ’  J 

W.W  s  aP^/A.^ 

i,l  1,1  i 


m 


(29 -a) 

(29-b) 


-1  1- 


Ww  „A[k]  /A  no 

i,m+l  i,m  1 


UV-c) 


The  ratio  of  measured  to  calculated  activity  is  defined  as 


R.W  -A./A.CkJ 
1  1  1 


Using  Eqs.  (29)  and  (30),  an  activity -weighted  correction  term  C,  (based 

rk]  rk  i  ^  . 

on  the  logarithm  of  R)'  since  the  R)  often  differ  by  several  orders  ol'  magni¬ 
tude  at  the  start  of  iteration,  depending  on  the  initial  approximation  spectrum) 

is  selected: 


Cj^  -  fn  R[k]/y^  jk^ 


for  all  j  such  that,£,  W.  ,  i  0.  Iteration  is  then  performed  by: 
1=1 


.  Ck+l]  _  ,*[k] 


HI 


For  all  j  such  that  £,  W,  <  =  0,  log-log  interpolation  is  performed  between 
i=l  1  *  J  n  fk"l 

the  nearest  lower  j  and  the  nearest  higher  j  for  which  £.  W,  ,  t  0  to  obtain  the 

1=1  1 J J 

next  (k+l)  value  of  differential  flux;  if  such  J  cannot  be  found  both  above  and 
below,  extrapolation  is  performed  (below  or  above,  as  required)  according  to 
one  of  several  alternate  forms: 

low  energy  end:  (1)  1/E; 

(2)  fission  (v^E)  fo  rm; 

(3)  Godiva  (E)  form, 

high  energy  end:  (1 )  fission  (Maxwellian  form) 

(2)  fusion  (14  Mev  peak  of  arbitrary  magnitude), 

t  h  , 

Iteration  is  repeated  until  a  solution  is  obtained  on  the  (  iteration,  Using 
Eqs,  (32)  and  re-substituting  the  (k  +  l)**1  result  into  the  (k  +  2)**1,  etc,  ,  the  f**1 
solution  may  be  written  as 
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r#i  w 

<pr  -  =  -  exp 


i-1  \ 

§‘* 


HO 


p  n 

where  <p.L0J  represents  the  initial  input  approximation, 

Tfl  r«i 

In  addition  to  providing  a  differential  flux  solution  <fil  J(E.)  m  <p,  J  at  621 

J  J 

energy  points,  the  SAND  II  code  provides  an  integral  flux  solution 


*-^(E  e  1 


asJ  yEj 
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(E)  dE 
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fei 

The  absolute  value  of  the  solution  is  obtained  by  renormalizing  <fl  (E)  so  that 
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where  Aj  L  is  the  normalized  activity  produced  by  unity  total  integral  flux, 
Calculated  values  of  saturated  activity  are  provided,  given  by 


r«i  r«i  f  r«i  3  *N  X  a  (E) 


dE 


36) 


Calculated  values  of  foil  reaction  energy  limits  of  sensitivity  and  E^  are  also 

provided,  defined  such  that  some  large  portion  a  (currently  95%)  of  the  activu- 
th 

tion  of  the  i  foil  is  produced  by  neutrons  of  energy  E  >  E,,  and  the  same  por¬ 
tion  a  1b  produced  by  neutrons  of  energy  E  <  E^,  Equation  (37)  implies  the 
definition  of  and  E^s 
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The  above  equations  are  solved  by  the  SAND  il  code,  which  has  been  written 
in  Fortran  IV  for  the  CDC  6600  computer  and  is  documented  in  Volume  II. 


SECTION  IV 

ANALYTICAL  STUDY  OF  THE  ITERATIVE  METHOD 


A,  Requirements  for  Solution  Achievement 
1,  General  Requirements 

The  unknown  measured  flux  spectrum  form  satisfies  (except  for  experi¬ 
mental  errors)  the  equations 

Aj  =  AjW  (i  =  1.2.  ,  .  .  n>  (38) 

where,  as  previously  defined,  the  A.'s  are  the  measured  saturated  activities. 

1  fk  1 

In  principle,  an  exact  solution  (Aj,  =  A.  for  all  "i")  can  be  obtained.  Such  a 
solution  will  not  be  unique,  however,  and  will  not  necessarily  be  more  useful 
than  an  approximate  solution  because  of  experimental  errors  in  A.  and  (j^(E), 

Gold  (Ref,  9)  has  stated  that,  in  general,  neither  an  exact  nor  even  an  approxi¬ 
mate  solution  is  required,  but  that  the  solution  should  properly  satisfy  the  set  of 
integral  equations  and  certain  necessary  conditions  imposed  by  physical  impli¬ 
cations;  he  calls  such  a  solution  an  appropriate  solution. 

Assuming  that  the  initial  approximation  (k  =  0)  is  not  itself  an  appropriate 
solution,  the  iterative  procedure  is  applied  until  Eqs.  (38)  are  satisfied,  to 
within  appropriate  limits  of  acceptable  error,  If  the  initial  approximation  spec¬ 
trum  is  properly  selected  and  if,  in  addition,  the  essential  physical  conditions 
have  been  properly  considered,  the  iterative  spectrum  for  which 

AjSA^  (39) 


for  all  i  will  be  an  appropriate  solution,  such  that 


*r,'(Ej.  1 1  »  »(Ej.  Em  , , 


(40) 


for  all  Ej  within  the  range  of  foil  measurements,  where 
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*(Ej, 


Em  +  1> 


and  «p(E)  is  the  true  (unknown)  spectral  form. 

2.  SAND  II  Requirements 

The  SAND  II  code  provides  a  "best  fit"  neutron  flux  spectrum  for  a  given 
input  set  of  infinitely  dilute  saturated  foil  activity  measurements,  without  inter¬ 
ruption  in  the  calculational  procedure  for  manual  calculation  or  decision  proc¬ 
esses.  SAND  11  accepts  as  input  information: 

a.  A  measured  set  of  infinitely  dilute  saturated  activities; 

b.  foil  covers  (cadmium,  boron  10  and  gold)  used,  if  any,  and  their 
thicknesses; 

c.  an  initial  spectral  approximation  (or  the  option  of  selection  from  the 
reference  spectrum  library); 

d.  solution  acceptability  criterion:  the  degree  to  which  all  621  differ¬ 
ential  flux  values  for  two  successive  iterations  must  agree  in  order 
to  consider  the  solution  appropriate; 

e.  discard  test  criterion:  the  number  of  standard  deviations  of  the 
measured  activity  of  a  given  foil  (calculated  from  all  foils  used)  from 
the  activity  calculated  from  the  solution  spectrum,  above  which  value 
the  foil  measurement  is  to  be  discarded  and  the  solution  revised  to 
account  for  its  exclusion. 

These  input  requirements  and  the  manner  in  which  they  affect  the  solu¬ 
tion  have  been  studied  for  both  analytical  and  experimental  cases,  The  results 
of  these  studies  will  be  considered  in  Section  IV-B  and  Section  V. 

B.  SAND  II  Results  for  Selected  Test  Cases 
1.  Summary  Discussion 

To  determine  the  accuracy  of  the  SAND  II  iterative  procedure,  test  cases 
were  studied  which  include  different  ratios  of  resonance  to  fission  to  fusion  neu¬ 
trons.  The  analytical  tests  consisted  of  attempts  to  ascertain  the  degree  to 
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which  SAND  II  is  capable  of  obtaining  integral  flux  and  "unfolding"  differential 
spectral  structure  under  both  ideal  and  non-ideal  conditions  —  that  is,  with  and 
without  errors  in  activity  and  cross  sections. 

For  different  kinds  of  test  spectra,  "exact"  activities  were  calculated 
using  the  current  cross  section  evaluations  and  attenuation  calculational  tech¬ 
niques,  Some  of  these  "true"  spectra  were  chosen  from  SAND  II  solutions  for 
the  tests  considered  in  Section  V,  while  others  were  selected  on  the  basis  of 
being  representative  of  different  types  of  environments.  The  "exact"  activities 
were  then  considered  to  be  "measured"  activities  and  used  (with  and  without 
added  errors)  as  input  to  SAND  II,  along  with  both  acceptable  and  unacceptable 
input  approximations,  to  obtain  iterative  solutions. 

Some  of  the  test  cases  considered  were  selected  by  Biggers  and  Mehl 
so  that  the  "true"  form  would  be  unknown  to  the  authors  until  after  a  suitable 
solution  had  been  obtained  (Ref.  13,  14),  This  technique  was  valuable  in  study¬ 
ing  the  effectiveness  of  the  SAND  II  code  to  use  the  reference  spectrum  library 
to  select  the  initial  approximation,  since  in  some  of  these  cases  no  information 
was  provided  as  to  the  type  of  environment,  (These  cases  provided  the  most 
severe  tests  of  the  "unfolding"  technique,  since  there  seldom  is  a  real  case  for 
which  no  information  is  available.) 

The  results  of  these  analytical  tests  were  in  general  excellent.  Given 
sufficient  foil  energy  coverage,  the  code  has  proved  capable  of  unfolding  the 
general  differential  envelope,  as  well  as  many  structural  details,  especially  for 
the  energy  range  above  0,1  Mev,  most  adequately  covered  by  threshold  defectors. 

The  test  cases  discussed,  in  whole  or  in  part  in  Sections  IV-B-2,  -3, 

>!< 

and  -4,  are  identified  as  follows: 

Test  1  —  Ground  Test  Reactor  (GTR)  core  spectrum  (reference  spectrum 
19,  Volume  IV) 

Test  2  —  Fusion  Type  Spectrum  (FTS) 

*Since  the  current  low-energy  detectors  in  the  SAND  II  code  do  not  have  differ¬ 
ent  cross  section  shapes  (all  are  essentially  1 /v)  in  the  thermal  region,  no 
attempt  was  made  to  determine  the  shape  of  the  thermal  spectrum  for  these 
tests.  The  SAND  II  code  simply  established  a  magnitude  for  the  input  iorm 
(usually  Maxwellian)  in  relation  to  the  higher  energy  component  of  the  spectrum, 
In  general,  therefore,  the  thermal  components  of  the  spectra  are  omitted  in  the 
discussions  of  both  the  analytical  and  experimental  results  that  follow. 
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Test  3  - 
Test  4  — 

Test  5  — 

Test  6  — 

Test  7  - 
Test  8  — 

Test  9  - 

Test  10  —  Solution  spectra  for  the  Epithermal  Critical  Experiment 
Laboratory  (ECEL)  fast  reactor  test. 

Test  11  —  Solution  spectrum  for  a  homogeneous  Water  Moderated 
Graphite  Reflected  (WMGR)  water  cooled  reactor  test. 

Test  12  —  Graphite  Moderated  and  Reflected  Water  cooled  (GMRW) 
reactor  test. 

Test  13  —  Time  of  Flight  Measurements  (TFM)  of  the  leakage  spectrum 
from  a  TRIGA  mockup  with  a  2"  water  reflector  (Ref.  15). 

2.  Appropriateness  of  Solution  and  Solution  Criterion 

For  a  flux  spectral  determination  covering  the  entire  energy  region,  the 
number  of  foil  reactions  used  will  in  most  cases  be  between  10  and  20.  There 
is,  therefore,  no  unique  exact  solution  to  Eqs,  (38),  and  no  unique  "best  fit" 
solution  iterative  limit.  Any  reasonable  iterative  procedure,  such  as  that  in¬ 
corporated  in  SAND  II,  may  be  used  that  fulfills  certain  obvious  requirements 
of  stability;  for  example,  it  must  not  yield  a  solution  which  diverges  to  infinity 
or  produces  negative  fluxes,  The  criterion  of  convergence  is,  however,  not  re¬ 
quired;  it  is  only  necessary  that  all  solutions  obtained  by  iterations  beyond  that 
which  yields  an  appropriate  solution  remain  within  a  fixed  neighborhood  of  that 
solution,  the  size  of  the  neighborhood  depending  on  acceptable  experimental 
error  limitations. 

The  first  solution  acceptability  criterion  studied  was  based  on  the  degree 
to  which  the  pertinent  foil  activities  calculated  from  the  solution  spectrum  must 
agree  with  the  input  measured  activities.  This  criterion  was  found  to  be 


‘i'These  were  the  numbers  originally  used  for  the  analysis  of  these  cases  by  a 
partially  automated  iterative  procedure  (see  Reference  3). 
tThe  Case  3  and  5  results  are  not  discussed  in  this  report. 


Biggers1  Case  1 
Riggers'  Case  ?/ 

Biggers'  Case  3 

:{( 

Biggers'  Case  4 
Biggers  and  Mehl's  Case  1 
Biggers  and  Mehl's  Case  2^ 
Biggers  and  Mehl's  Case  4^ 


(Forms  assumed  to  be  unknown  to  authors, 
i,  e.  ,  used  SAND  II  code  with  its  reference 
spectrum  library  as  a  "black  box.") 


Forms  initially  unknown  to 
authors 
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unacceptable  for  most  test  cases  because  of  its  obvious  inherent  lack  of  sensi¬ 
tivity  to  differential  flux  changes,  A  series  of  tests  was  initiated,  therefore,  to 
study  this  parameter  as  well  as  other  aspects  of  the  iterative  procedure. 

A  parametric  study  of  the  formula  used  to  calculate  the  weighting  function, 
Eqs,  (29),  was  conducted  in  an  effort  to  minimize  what  at  first  appeared  to  be  in¬ 
appropriate  perturbations  in  the  solution  differential  flux  spectrum  over  certain 
energy  regions,  in  ill-conditioned  cases  which  include  large  errors  in  measured 
foil  activity  and/or  cross  sections. 

Results  of  this  study  made  it  clear  that  the  difficulty  (which  is  usually 
negligible  and  is  in  any  case  not  serious  from  the  standpoint  of  integral  flux  so¬ 
lution  accuracy)  is  not  one  of  method.  The  difficulty  is,  rather,  inherent  in  the 
problems  of  foil  coverage  and  activity  measurement,  non-uniqueness  of  solution, 
and  incomplete  or  inaccurate  evaluations  of  reaction  cross  section  data.  Physi¬ 
cally  unreasonable  structure  in  solution  spectra  can,  therefore,  be  minimized 
by  careful  experimental  planning  and  execution  to  assure  sufficient  foil  energy 
coverage  and  accuracy  of  activity  measurements,  by  utilization  of  all  available 
physical  information  to  select  an  appropriate  spectrum,  and  by  the  acquisition 
of  sufficiently  accurate  cross  section  data  (see  the  results  for  Tests  10,  11, 
and  1  2  in  Section  V-B). 

The  appearance  of  inappropriate  perturbations  when  inaccurate  cross 
section  data  are  used  is  actually  an  important  advantage,  since  it  can  be  used 
to  survey  currently  used  cross  section  evaluations  to  determine  which  require 
revision  and  in  what  energy  regions  the  revision  is  needed.  The  SAND  II  code 
could  then  be  adapted  to  provide  improved  cross  section  evaluations  as  solutions. 
Because  of  the  "best  fit"  approach  of  SAND  II  which  provides  a  solution  using 
redundant  (in  energy  coverage)  foil  reactions,  the  solution  behavior  can  also  be 
examined  to  identify  inordinately  large  activity  measurement  errors,  thus  pro¬ 
viding  a  tool  to  help  identify  either  systematic  errors  in  determining  activities 
or  uncorrected  absorption-scattering  effects  (see  the  results  for  Test  10  in 
Section  V - B), 

A  "recycle"  approach  was  utilized  to  help  in  the  study  of  the  appropriate¬ 
ness  of  a  given  iterative  solution.  This  was  done  by  considering  the  SAND  II 
calculated  values  of  activity  for  a  solution  spectrum  to  be  a  new  set  of  exact 


"measured"  activities.  The  solution  spectrum  was  then  used  in  the  recycle  test 
as  the  "true"  spectrum  with  the  new  set  of  exact  "measured"  activities  and  the 
same  input  spectral  approximation  (that  in  each  case  had  led  to  the  "true"  spec¬ 
trum)  as  input  to  SAND  II. 

This  recycle  procedure  helps  to  identify  those  features  of  a  solution 
spectrum  which  are  caused  by  peculiarities  of  the  iteration  technique  or  any 
other  aspect  of  the  data  reduction  method;  if  no  such  peculiarities  exist  to  cause 
erroneous  spectral  structure,  the  solution  to  a  test  using  this  procedure  should 
agree  with  the  "true"  spectrum.  (It  is  recommended  that  this  procedure  be 
applied  to  any  real  case  in  which  the  reality  of  spectral  structure  comes  into 
question.)  A  typical  test  result  is  shown  in  Figure  1  for  Test  12.  The  100- 
iteration  limit  was  chosen  to  test  whether  once  a  reasonable  solution  had  been 
achieved,  further  iteration  would  cause  the  solution  to  diverge,  or  "blow  up." 

In  all  cases  studied,  it  was  found  that,  once  an  appropriate  solution  is  achieved, 
it  remains  within  a  quite  acceptable  neighborhood  of  that  solution  upon  further 
iteration. 

It  was  found  that  certain  ill -conditioning  can  cause  an  apparent  "blow  up" 
(depending  on  the  method  of  presentation  of  results).  This  occurs  when  a  reac¬ 
tion  has  a  response  curve  such  that  it  is  very  dominant  in  a  particular  narrow 
energy  region,  compared  to  all  other  reactions  used,  and  the  measured  activity 
for  that  reaction  is  in  substantial  error  (again  compared  to  the  other  reactions) 
on  the  low  side.  This  causes  negative  corrections  to  the  flux  at  that  energy  re¬ 
gion;  if  the  first  few  iterations  are  not  sufficient  to  generate  close  agreement 
between  the  measured  and  calculated  activities  for  that  reaction,  subsequent 
iterations  will  cause  continued  negative  corrections  to  the  flux  at  that  energy 
region.  After  the  first  several  reductions,  however,  which  have  been  sufficient 
to  remove  most  of  the  flux  from  that  region,  the  continued  negative  corrections 
although  perhaps  substantial  relatively,  will  have  no  further  significant  correc¬ 
tive  effect  on  the  calculated  activity  (which  is  produced  over  a  much  wider  energy 
range,  not  so  strongly  affected  by  the  negative  error  of  that  foil).  If  the  spectrum 
is  graphically  plotted  on  a  logarithmic  display  (as  is  usual),  successive  iterations 
will  cause  a  successively  more  severe  "dip"  at  the  pertinent  energy  region  while, 
actually,  once  the  flux  is  "almost"  zero,  compared  to  the  flux  in  its  immediate 
energy  vicinity,  further  reduction  (which  appears  substantial  on  a  logarithmic 
display)  will  have  no  significant  effect. 
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NEUTRON  DIFFERENTIAL  FLUX  (n  cnK-sec  Mev! 
(ARBITRARY  SCALE) 


It  is  emphasized  that  this  apparent  (on  a  logarithmic  display)  blow  up 
can  occur  only  in  fairly  ill-conditioned  cases,  when  a  set  of  foils  is  used  such 
that  one  foil  (1)  is  particularly  dominant  in  one  narrow  energy  region,  (2)  lias 
an  activity  produced  over  a  much  wider  energy  range,  over  most  of  which  it  is 
not  particularly  dominant  compared  to  the  other  foils  used,  and  (3)  has  a  large 
negative  measured  activity  error. 

In  the  resonance  region,  the  code  cannot  always  produce  real  structure 
in  the  solution  spectrum  when  the  resonance  foils  are  used  in  the  usual  manner; 
conversely,  the  code  will  often  produce  artificial  structure.  This  is  because 
the  very  severe  resonance  structure  in  the  reaction  cross  section  functions  can 
result  in  sensitivity  energy  ranges  for  the  various  resonance  foils  which  are  es¬ 
sentially  separate.  In  some  extreme  cases  (notably  gold),  90%  or  more  of  the 
activity  can  be  produced  by  one  major  cross  section  resonance  peak.  This 
means  that  initial  iterative  flux  corrections  are  generally  not  averaged  over 
more  than  one  foil,  and  peaks  and  dips  are  readily  introduced  at  the  respective 
resonance  energies,  to  make  the  spectrum  consistent  with  all  the  measured  ac¬ 
tivities.  This  difficulty  will  not  be  encountered  in  a  spectrum  with  very  little 
thermal  and  intermediate  flux,  since  this  would  result  in  response  curves  shifted 
upward  in  energy  from  the  resonance  region  into  a  region  in  which  they  are  much 
smoother  and  overlap,  Usually,  only  1/E  or  similar  behavior  in  the  intermediate 
energy  range  will  allow  the  problem  to  arise.  As  previously  stated,  the  diffi¬ 
culty  of  artificial  structure  may  in  some  cases  be  put  to  good  advantage;  exami¬ 
nation  of  the  solution  lot  structure  of  this  type  is  often  an  easy  way  to  detect 
large  errors  in  activity  measurements  or  self-absorption  effects  (see  Sec¬ 
tion  V-B). 

A  very  effective  device  for  solution  of  the  problem  of  artificial  structure 
in  the  resonance  region  is  the  use  of  self  -  sandwiched  resonance  detectors.  This 
would  preferentially  attenuate  a  large  portion  of  the  fiux  at  the  resonance  ener¬ 
gies  in  the  outer  covers,  effectively  smoothing  the  response  curves  for  the  sand¬ 
wiched  foils  by  eliminating  or  sharply  reducing  the  large  peaks;  this  in  turn 
causes  regions  of  foil  energy  overlap  to  be  much  more  significant  in  terms  of 
contribution  to  total  activity.  In  effect,  self-sandwiching  causes  resonance  foils 
to  behave  more  like  integrating  threshold  detectors. 
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The  results  of  he  above  studies  as  well  as  those  that  wili  be  discussed 
in  subsequent  sections  of  this  report  suggest  that  the  criterion  for  solution 
achievement  should  be  based  on  a  comparison  of  the  differential  spectra  of  suc¬ 
cessive  iterations.  For  the  current  SAND  II  code,  therefore,  a  particular  itera¬ 
tive  spectrum  is  considered  to  be  an  appropriate  solution  if  its  differential  flux 
differs  from  that  of  the  preceding  iteration  by  less  than  a  given  percent  (speci¬ 
fied  as  input)  at  all  621  energy  points  between  1x10  ^  Mev  and  18  Mev. 

The  results  of  a  parametric  study  of  the  required  number  of  iterations 
as  a  function  of  the  value  of  the  current  solution  criterion  (as  well  as  of  the  study 
based  on  the  standard  percent  deviation  of  the  ratios  of  measured  to  calculated 
activity),  for  a  measurement  of  the  neutron  flux  spectrum  in  the  ECEL  reactor, 
are  presented  in  Figure  2.  (Other  aspects  of  the  ECEL  test  are  discussed  in 
detail  in  Section  V-B,  )  Testing  performed  to  date  suggests  that,  in  most  realis¬ 
tic  cases,  a  successive  iteration  difference  of  about  5  to  10%  is  a  reasonable 
solution  criterion,  and  that  such  a  criterion  is  usually  satisfied  in  fewer  than 
thirty  iterations.  The  actual  value  used,  however,  will  depend  on  the  confidence 
the  user  places  in  the  accuracy  of  the  measured  activities,  and  on  the  aptness  of 
the  initial  spectral  approximation. 

3,  Foil  Set  Selection  and  Input  Spectrum  Sensitivity 

During  the  development  and  testing  of  the  SAND  II  code,  studies  were 
performed  to  determine  the  degree  of  input  sensitivity  of  the  current  iterative 
procedure;  that  is,  to  what  degree  the  solution  differential  flux  shape  is  affected 
by  the  selection  of  the  foil  set  and  the  choice  of  the  input  spectral  approximation, 
Based  on  the  studies  completed,  it  appears  that,  if  the  foil  set  selected  provides 
sufficient  coverage,  the  code  is  less  input  dependent  than  was  originally  expected 
In  cases  where  foil  coverage  is  not  adequate,  however,  the  code  has  been  ob¬ 
served  to  be  quite  input  sensitive  in  terms  of  the  differential  flux.  A  poor  choice 
of  input  spectrum  can  result  in  foil  sensitivity  energy  regions  which,  if  insuffici¬ 
ent  foil  overlap  occurs,  cannot  be  shifted  sufficiently  by  iteration  to  approach 
the  true  solution,  As  suggested  previously,  a  "recycle"  technique  can  be  used  to 
study  the  input  sensitivity  of  the  iterative  procedure,  in  those  cases  in  which  the 
credibility  of  the  solution  spectrum  is  in  question. 
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The  dependence  of  the  SAND  II  solution  spectrum  on  the  selection  of  a 

suitable  foil  set  and  an  initial  spectral  approximation  can  beat  be  illustrated  by 

presenting  several  examples,  The  SAND  II  solution  spectrum  will  not  agree 

with  the  "true"  spectrum  below  the  energy  range  of  the  lowest  detector  (about 

10  Mev  for  most  of  the  tests  except  Test  7,  for  which  it  is  about  10  Mev 

-4 

because  of  the  sharp  spectral  cutoff  at  about  10  Mev);  the  complete  spectrum 
is  retained  in  Tests  7,  8,  and  9  simply  for  completeness. 

First,  several  cases  were  studied  in  which  a  large  number  (~-25  of  the 
29  available  reactions)  of  foils  were  used,  and  no  information  was  provided  on 
the  errors  associated  with  the  measurement  of  the  foil  activities  and/or  the  type 
of  neutron  environment  being  measured.  The  only  available  information  was 
that,  in  two  of  five  cases  (three  of  which  are  discussed  here  as  Tests  7,  8, 
and  9),  foil  activity  probable  errors  are  10%;  in  one  case  (Test  7),  an  initial  in¬ 
put  spectral  form  was  provided  (shown  by  the  dashed  curve  in  Figure  3), 

The  SAND  I.T.  results  for  these  three  cases,  and  the  "true"  spectral  forms, 
are  given  in  Figures  3,  4,  and  5  respectively.  The  reference  spectrum  library 
was  used  to  select  an  initial  approximation  for  Tests  8  and  9.  For  Test  9,  the 
library  search  technique  was  modified  to  include  two  separate  reference  library 
searches,  one  for  the  low  energy  (resonance)  and  one  for  the  high  energy  (thresh¬ 
old)  foil  detectors.  The  two  spectra  selected  were  then  combined  to  provide  a 
hybrid  input  approximation.  (The  application  of  this  technique  can  significantly 
increase  the  confidence  ir.  the  final  SAND  II  solution  when  the  reference  spec¬ 
trum  library  is  used,  since  it  results  in  an  initial  approximation  with  a  reduced 
standard  deviation  of  activity  ratios,)  The  mutually  independent,  selection  of  two 
regions  of  the  initial  approximation  spectrum  is  justified  for  spectra  in  which 
the  resonance  foils  and  threshold  foils  form  essentially  independent  sets  —  that 
is,  their  sensitivity  curves  do  not  substantially  overlap  in  energy. 

In  general,  the  differential  flux  results  are  reasonably  correct  for  these 
three  tests,  even  though  the  foil  coverage  (and  overlap)  was  rather  poor  in  cer¬ 
tain  energy  regions.  This  lack  of  proper  foil  coverage  might  seem  surprising 
in  view  of  the  large  number  of  foils  used  for  these  tests.  It  must  be  emphasized, 
however,  that  the  SAND  II  energy  range  spans  more  than  eleven  decades;  also, 
the  25  foil  reactions  used  do  not  necessarily  have  energy  response  functions 
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NOTE:  THE  SAND  II  CODE  SELECTED  REFERENCE  SPECTRUM  29  (VOLUME  IV) 
AS  THE  INITIAL  INPUT  APPROXIMATION. 
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Figure  4.  SAND  II  Differential  Fluen.ce  Result  for  Test  8 


NOTE:  FOR  THIS  TEST  THE  LOW-  AND  HIGH-ENERGY  DETECTORS  WERE 
SEPARATED  TO  ALLOW  THE  SAND  II  CODE  TO  MAKE  BOTH  A  LOW- 
AND  HIGH-ENERGY  SEARCH  WITH  THE  REFERENCE  LIBRARY;  AS 
A  RESULT  OF  THESE  SEPARATE  SEARCHES,  AN  INPUT  APPROXI¬ 
MATION  WAS  SELECTED  THAT  IS  A  HYBRID  OF  REFERENCE 
SPECTRA  51  AND  23  (VOLUME  IV), 
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Figure  5.  SAND  II  Differential  Fluence  Result  for  Test  9 


suitable  for  obtaining  a  proper  solution  in  all  cases  especially  if  selected  with 
little  or  no  knowledge  of  the  environment. 

SAND  II  integral  flux  results  for  similar  cases  (Tests  3,  4,  5,  and  6) 
will  be  discussed  in  more  detail  in  Section  IV-B-4  for  a  smaller  set  of  14  foil 
detectors  that  have  been  found  suitable  for  measurements  in  a  number  of  fission- 
and  fusion-type  environments. 

The  next  pertinent  examples  consider  only  the  high  energy  region  of  the 
neutron  spectrum,  above  ~0.1  Mev.  The  iterative  results  of  two  "recycle" 
studies,  Tests  11  and  13,  are  presented  in  Figure  6.  The  same  set  of  13  foil 
detectors  was  used  for  Test  13A  and  Test  11.  Test  1  3B  used  the  same  set  of 
15  foils  that  was  used  for  Test  12,  Figure  1,  The  13-foil  set  included  sulfur, 
while  the  15-foil  set  did  not.  (The  individual  foils  used  for  these  tests  are  listed 
in  Tables  XIV  and  XVIII  of  Section  V.) 

In  Figure  6,  the  "true"  spectrum  and  SAND  II  results  for  each  test  are 
represented  by  the  solid  lines  and  the  input  approximation  by  the  dashed  line, 
Tests  11  and  13  again  demonstrate  that  the  iterative  method  is  capable  of  unfold¬ 
ing  differential  structure,  but  both  the  relative  hardness  of  the  spectrum  (which 
will  affect  the  energy  response  range  of  individual  foils),  and  the  particular  foil 
set  selected,  will  affect  the  final  shape  of  the  solution  differential  flux  spectrum, 
(The  effect  on  the  integral  flux  is  significantly  less.)  In  this  regard,  it  is  noted 
that  the  general  envelope  of  the  solution  differential  flux  is  more  correct  in 
Test  12,  Figure  1,  than  in  either  Tests  11  or  13.  This  is  at  least  in  part  due  to 
the  fact  that  in  Test  12  the  solution  spectrum  is  softer  and  more  (fifteen)  foil 
detectors  were  used.  The  fine  structure  superimposed  on  the  solution  envelope 
between  2.5  and  5  Mev  for  Test  11  is  due  mainly  to  similar  fine  structure  in  the 
cross  section  of  the  S  (n,p)P  reaction  (see  Volume  III),  It  is  most  pro¬ 
nounced  in  this  test  because  of  the  relative  softness  of  the  spectrum  and  the 
particular  foil  set  used.  Conversely,  the  effect  of  sulfur  is  only  slightly  evident 
in  Test  13A  where  the  spectrum  is  much  harder.  The  sulfur  cross  section  fine 
structure  is  predominantly  in  the  low  energy  region  of  its  sensitivity  range; 
therefore,  a  softer  environment  causing  increased  low-region  sensitivity  and 
therefore  more  dominant  low-end  iterative  weighting  (see  Eqs.  (29)).  will  result 
in  a  stronger  reflection  of  this  structure  in  the  solution  spectrum  than  will  a 
harder  environment. 
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Figure  6,  SAND  II  Recycle  Study  Differential  Flux  Results 
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NEUTRON  DIFFERENTIAL  FLUX  (n.W-sec-Mev) 
(ARBITRARY  SCALE) 


These  examples  serve  to  illustrate  the  dependence  of  the  input  sensitivity 
of  the  iterative  solution  on  the  choice  of  foil  sets  and  initial  spectral  approxima¬ 
tion,  A  final  example  illustrates  a  case  in  which  all  the  high  energy  foils  were 
utilized  and  three  different  input  spectral  approximations  were  selected.  Fig¬ 
ure  7  presents  the  iterative  results  for  Test  2.  An  assumed  fusion  form  was 
selected  as  the  "true"  spectrum  for  this  test;  the  three  different  input  approxi¬ 
mations  (the  three  dashed  lines)  were  used  to  obtain  the  three  solutions  (three 
solid  lines)  shown.  This  test  makes  apparent  the  capability  of  the  SAND  II  code 
to  generate  a  reasonable  solution  with  a  poor  input  approximation  in  certain  cases 
when  there  is  adequate  foil  coverage, 

4.  Integral  Flux  Comparison 

As  indicated  previously,  the  thermal  component  of  the  flux  has  been 
omitted  in  the  discussion  of  differential  flux  results  throughout  this  report,  be¬ 
cause  currently  available  low  energy  detectors  do  not  have  sufficiently  different 
response  functions  to  allow  spectral  unfolding.  The  thermal  component  of  the 
integral  flux  spectra  will  not  be  considered  either,  since  all  that  is  currently 
possible  is  to  establish  the  magnitude  of  the  flux  in  this  region. 

It  is  worth  mentioning,  however,  that  there  are  foil  reactions  that  have 
only  one  predominant  resonance  in  the  few  electron  volt  region.  Combinations 
of  such  foils  (both  bare  and  cadmium-covered)  present  a  variety  of  absorption 
cross  sections,  which  could  be  used  to  determine  the  shape  of  the  differential 
flux  in  the  thermal  region.  Fulmer  and  Ruane  (Ref,  16)  have  recently  reported 
on  the  use  of  such  a  group  of  detectors  for  determining  spectra  in  this  region. 

The  incorporation  of  such  reactions  as  Dy*^(n,y  )Dy*^,  Eu*  ^*(n,y)Eu* 
Lu^^(n,y)Lu^^,  and  In^  ^(n,y)In*^  into  the  SAND  II  cross  section  library 
could  permit  the  unfolding  of  spectra  for  the  very  low  energy  region. 

To  determine  the  accuracy  of  the  SAND  II  integral  flux  results,  several 
test  cases  were  studied.  Test  1  involved  a  known  spectrum  for  the  Ground  Test 
Reactor  (reference  spectrum  19),  a  set  of  foil  activities  with  and  without  errors, 
and  the  use  of  a  previously  determined  input  approximation  spectrum,  (This  was 
a  repeat  of  a  test  that  had  been  studied  before  and  was  reported  in  reference  2.) 
The  SAND  II  results  are  presented  in  Table  I  and  Figure  8  for  the  first  part  of 
this  test,  which  involved  the  assumption  that  there  were  no  errors  in  measured 
activities  or  cross  sections. 


-32- 


NOTE:  THE  TRUE  SPECTRUM  IS  ARBITRARILY 
NORMALIZED  AT  1.0  Mev 
FOR  EACH  TEST  SHOWN. 


INPUT  SPECTRUM:  CONSTANT. 
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SPECTRUM  31  (VOLUME  IV). 
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Figure  7.  SAND  II  Differential  Fluence  Result  with  Different 
Input  Approximations  for  Test  2 
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NEUTRON  DIFFERENTIAL  FLUENCE  (r.  raf-Mev') 
(ARBITRARY  SCALE) 


ENERGY  (Mev) 

D  II  Iterative  and  True  Test  Spectrum  Integral  Flux  Comparison 
Ground  Test  Reactor  (GTR)  Core  Flux  Spectrum  Test  1  with 


The  results  for  Test  1  were  intended  to  be  representative  of  a  typical 
well- moderated  reactor  test  where  the  user  of  the  SAND  II  code  might  have 
some  pre  -  knowledge  of  the  anticipated  overall  form  of  the  unknown  spectrum. 

It  was  assumed,  therefore,  as  might  be  in  an  actual  case,  that  a  1  /  E  +  fission 
spectrum  (reference  library  spectrum  7)  is  a  reasonable  input  approximation 
based  on  the  available  physical  knowledge  about  the  reactor  environment.  'More 
typically,  a  good  input  approximation  might  be  available  as  a  result  of  reactor 
physics  calculations.) 

The  results  for  Test  1  are  good,  considering  that  the  foil  set  selected 
did  not  provide  proper  foil  coverage  for  the.  regions  ~10  ^  to  1  Mev  and  1  5 
to  18  Mev  (see  Table  I),  where  negative  accumulated  errors  up  to  and  greater 
than  •''•19%  appear  in  Figure  8.  These  large  negative  errors  are  a  result  of  a 
lack  of  proper  foil  coverage,  which  yields  a  solution  that  is  largely  a  reflection 
of  the  input  approximation  in  these  regions.  For  those  regions  where  there  is 
adequate  foil  coverage,  however,  the  integral  flux  solution  is  accurate  to  within 
±10%  based  on  a  relative  comparison;  i.  e.  ,  neglecting  the  "-19%  error  accumu¬ 
lated  between  10  *  and  1  Mev.  Similar  results  are  given  in  Table  II  and  Fig¬ 
ure  9  for  the  second  part  of  this  test,  where  activity  errors  were  introduced, 
selected  at  random  from  a  set  of  percent  error  values  with  an  assumed  normal 
distribution  and  a  standard  deviation  of  15%.  In  Figure  9,  the  manner  .in  which 
the  integral  flux  solution  reflects  the?  magnitude  and  the  direction  of  the  errors 
listed  in  Table  II  is  apparent.  In  this  case,  the  integral  flux  solution  is  within 
±20%  at  most  energies.  (Experimental  results  for  environments  similar  to 
that  of  Test  1,  where  the  initial  input  spectra  were  based  on  reactor  physics 
calculations,  are  discussed  in  Section  V-B.) 

A  cooperative  study  was  conducted  involving  work  performed  at  Atomics 
International  and  the  Los  Alamos  Scientific  Laboratory.  For  the  initial  stop  of 
the  study,  the  evaluated  cross  section  data  presented  in  reference  17  were  used 
by  Diggers  of  Los  Alamos  Scientific  Laboratory  to  calculate,  for  each  of  four 
selected  spectra  .  the  activities  of  a  set  of  14  foils  selected  from  the  29  reactions 
listed  in  Volume  III.  The  calculated  saturated  activities  of  the  selected  foils  for 
each  of  the  four  test  cases  wt  re  then  altered  in  a  random  manner  to  simulate 
experimental  errors,  and  were  used  as  "measured"  activities  for  the  test  cases. 
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Random  Activity 


The  selection  of  a  best  set  of  foil  detectors  depends  on  having  adequate 
knowledge  about  the  anticipated  form  of  an  unknown  spectrum.  Without  such 
knowledge,  a  selected  set  may  not  contain  foils  with  suitable  energy  sensitivity 
limits  for  the  unknown  spectrum.  Since  the  forms  of  the  4  test  spectra  were 
originally  not  known,  the  set  of  14  foil  reactions  used  was  based  on  a  set  that 
had  previously  been  found  suitable  for  both  fission  and  fusion  type  environments. 

These  four  cases  were  used  to  test  the  capability  of  the  SAND  II  code  to 
select  a  proper  input  approximation  from  the  current  reference  library.  The 
results  presented,  therefore,  are  representative  of  the  type  of  solution  that 
might  be  obtained  when  the  SAND  II  code  is  used  in  a  "black  box"  fashion. 

Biggers  had  selected  four  spectra  with  different  ratios  of  resonance  to 
fast  to  fusion  neutrons,  identified  as  Tests  3,  4,  5,  and  6.  He  had  used  the 
cross  section  data  of  reference  17  to  calculate  the  saturated  activities  of  the  14 
foils,  for  each  of  the  four  test  cases.  These  calculations  utilized  a  group  struc¬ 
ture  with  four  energy  groups  per  decade  (group  energy  limits  within  each  decade: 
1,  2,  4,  7)  below  1.0  Mev,  and  0.5  Mev  intervals  above  1.0  Mev,  Group- 

>!< 

averaged  cross  sections  based  on  uniform  flux  weighting  within  each  group  were 
used  to  calculate  the  contribution  of  each  group  to  the  saturated  activity  of  each 
foil.  The  results  were  summed  for  each  foil  and  a  random  error  was  introduced 
into  each  of  the  14  saturated  activity  values.  The  activities  and  assigned  values 
of  error  as  provided  by  Biggers  for  the  4  test  cases  are  given  in  Table  III, 

For  these  tests,  the  solution  criterion  and  foil  discard  values  were  as¬ 
signed  extreme  values  so  that  no  foils  would  be  discarded  and  the  iteration  would 
proceed  to  an  assigned  limit  of  30  cycles.  Although  this  procedure  is  unrealistic 
from  an  experimental  standpoint,  it  does  make  the  effect  of  experimental  errors 
(those  assigned  in  Table  III)  more  pronounced  for  studying  the  iterative  process. 

In  a  completely  "black  box"  fashion,  the  SAND  II  code  selected  the  initial 
spectral  approximation,  performed  the  numerical  steps  of  iteration,  and  pro¬ 
vided  a  solution.  The  detector  results  for  Tests  3,  4,  5,  and  6  are  presented 
in  Tables  IV,  V,  VI  and  VII  respectively.  An  integral  flux  comparison  was  made 
of  the  SAND  II  solutions  with  the  "true"  spectra  for  all  cases,  and  the  results  are 

"This  is  a  substantially  more  coarse  approximation  than  that  used  in  SAND  II, 
where  the  energy  interval  division  is  much  finer  (see  Section  III-A), 
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Table  IV 

SAND  II  SOLUTION  RESULTS  OBTAINED  AFTER  30  ITERATIONS  FOR  TEST  3 


41 


Standard  Deviation  of  Measured  Activities  4.92% 


Table  VI 

SAND  n  SOLUTION  RESULTS  OBTAINED  AFTER  30  ITERATIONS  FOR  TEST  5 
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*Means  foil  detector  is  covered  with  cadmium. 


SAND  II  SOLUTION  RESULTS  OBTAINED  AFTER  30  ITERATIONS  FOR  TEST  6 


given  in  Tables  VIII,  IX,  X,  and  XI.  The  results  for  Test  3  are  also  displayed 
graphically  in  Figure  10  for  comparison  and  further  discussion. 

For  Test  3,  and  as  explained  in  reference  3,  no  solution  could  be  ob¬ 
tained  with  the  SAND  I  code  "by-hand"  procedure  below  1  Mev,  because  the  foil 
set  selected  was  not  suitable  for  this  spectrum.  With  SAND  II,  however,  a 
satisfactory  solution  was  obtained  and,  in  fact,  the  value  of  total  flux  obtained 

is  seen  to  be  very  accurate  (less  than  1/4%  error).  The  relatively  large  error 

.  3 

(nearly  20%)  in  the  region  about  10  Mev  probably  results  from  an  assumed 
fission  form  in  the  low  energy  range  for  the  input  approximation,  while  the  true 
spectrum  had  zero  flux  below  10"^  Mev  (the  arbitrary  low  limit  of  the  calcula¬ 
tion  by  Biggers).  The  "blowup"  of  the  solution  in  the  region  above  approximately 
14  Mev  results  from  the  same  kind  of  pathology  in  the  true  spectrum,  which  was 
assigned  zero  flux  above  15.5  Mev.  This  type  of  behavior  was  not  found  in  other 
tests,  where  there  were  no  assigned  sharp  flux  cutoff  points.  No  provision  is 
included  in  SAND  II  for  coping  with  spectra  having  a  sudden  cutoff  limit. 

The  effect  of  experimental  errors  on  the  solution  was  amplified  in  these 

tests  by  setting  the  solution  criterion  value  very  low,  <0,1%,  so  that  iteration 

would  continue  to  the  specified  limit  of  30  cycles;  this  allows  the  integral  solu- 

+ 

tion  to  develop  more  structure  to  become  more  nearly  exact.  The  effect  of 
these  errors  on  the  solution  is  clearly  seen  in  a  comparison  of  the  results  in 
Table  XI  for  the  energy  region  between  4  and  8  Mev  with  the  results  for  the  in- 
dividual  foils,  Tables  III  and  VII,  for  Test  6,  The  A1  (n,  p)  assigned  value  of 
error  (-13.1%)  is  obviously  in  error  itself.  The  deviation  of  measured  from 
calculated  activity  for  this  one  reaction  is  -59.98%,  which  is  2.42  standard  devi¬ 
ations,  Based  on  a  reasonable  foil  discard  criterion  of  1,98  standard  deviations 
(95%  confidence  level  assuming  normally  distributed  errors),  this  foil  would 
have  been  discarded  in  a  normal  SAND  II  run  (see  Section  IV-A-2).  In  a 


>!<The  SAND  I  iterative  procedure  involved  graphical  plotting  of  integral  flux 
points  and  differentiation  by  a  "by-hand"  best  fit  integral  flux  curve. 
fExactness  in  this  sense  refers  to  the  degree  of  consistency  of  the  input 
"measured"  activities,  or  the  degree  to  which  the  ratios  of  measured  to  cal¬ 
culated  activities  for  the  spectrum  in  question  agree;  exactness  is  not  neces¬ 
sarily  meant  to  infer  correctness. 
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Table  VIII 


TEST  3  -  COMPARISON  OF  THE  SAND  II  30th  ITERATION 
ABSOLUTE  INTEGRAL  FLUX  RESULTS  WITH 
THE  TRUE  TEST  SPECTRUM 


>!< 

Neutron 

Energy 

(Mev) 

Integral.  Neutron  Flux 

Ratio 

SAND  II /TTS 

True  Test 
Spectrum 
(TTS) 

SAND  II 
Spectrum 

4  x  10"7 

- 

3.171  x  104 

-  4 

1  x  10 

.  + 

3.169  x  104 

- 

1  x  10”3 

3.282  x  104 

2.840  x  104 

0.865 

1  x  10”2 

2.167  x  104 

2.300  x  104 

1.061 

1  x  10”1 

1.786  x  104 

1.509  x.  104 

0.845 

1.0 

9.343  x  103 

9.354  x  103 

1.001 

1.5 

6.963  x  103 

7.032  x  103 

1.010 

2.0 

5.063  x  103 

5.155  x  103 

1.018 

3.0 

3.055  x  103 

3.138  x  103 

1.027 

4.0 

2.172  x  103 

2.176  x  103 

1.002 

5.0 

1.578  x  103 

1.676  x  103 

1.062 

6.0 

1.280  x  103 

1.349  x  103 

1.054 

7.0 

1.112  x  103 

1.122  x  103 

1.009 

8.0 

9.935  x  102 

9.228  x  102 

0.920 

9.0 

9.122  x  102 

7.866  x  102 

0.862 

10.0 

8.306  x  102 

7.158  x  102 

0.862 

11.0 

7.548  x  102 

6.804  x  102 

0,901 

12.0 

6.589  x  102 

6.504  x  102 

0.986 

13.0 

5.273  x  102 

5.927  x  102 

1.124 

14.0 

2.503  x  102 

3.063  x  102 

1.224 

15.5 

0 

1.971  x  10 1 

- 

’i'Not  all  of  the  "true"  spectrum's  energy  points  are  li  ted.  j 
tThe  "true"  differential  flux  was  set  equal  to  zero  below  10”'  Mev. 
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Table  IX 


TEST  4  -  COMPARISON  OF  THE  SAND  !f  30th  ITERATION 
ABSOLUTE  INTEGRAL  FLUX  RESULTS  WITH 
THE  TRUE  TEST  SPECTRUM 


Integral  Neutron  Flux 

Ratio 

SAND  11/ TTS 

True  Test 
Spectrum 
(TTS) 

SAND  II 
Spectrum 

4  x  10-7 

1.325  x  105 

_ 

1  x  10-4 

1 

tit  ' 

1.325  x  105 

- 

1  x  10"3 

1.570  x  10 

1.322  x  105 

0.842 

4  x  10‘3 

1.565  x  105 

1.316  x  105 

0.841 

1  x  10"2 

1.555  >:  105 

1.306  x  105 

0.840 

4  x  10"2 

1.541  x  105 

1.294  x  105 

0.840 

1  x  10-1 

1.521  x  105 

1.289  x  105 

0.847 

4  x  10” 1 

1.431  x  105 

1.284  x  105 

0.897 

1.0 

1.251  x  105 

1.218  x  105 

0.947 

1.5 

1.201  x  105 

1.121  x  105 

0.933 

2.0 

1.126  x  105 

1.039  x  105 

0.921 

3.0 

8.920  x  104 

9.036  x  104 

1.013 

4.0 

8.525  x  104 

7.902  x  104 

0.927 

5.0 

7.550  x  104 

7.000  x  104 

0.927 

6.0 

6.600  x  ,104 

6.300  x  104 

0.955 

7.0 

5.830  x  104 

5.791  x  104 

0.993 

8.0 

5.430  x  104 

5.369  x  104 

0.989 

9.0 

4.755  x  104 

4.983  x  104 

1.048 

10.0 

4.155  x  104 

4.626  x  104 

1.114 

11.0 

3.775  x  104 

4.259  x  104 

1.128 

12.0 

3.435  x  104 

3.864  x  1()4 

1.125 

13.0 

3.135  x  104 

3.275  x  104 

1.045 

14.0 

2.020  x  104 

2,009  x  104 

0.995 

15.0 

2.900  x  103 

2.494  x  103 

0.860 

>!<Not  all  of  the  true  spectrum's  energy  points  are  listed,  ^ 
fThe  "true"  differential  flux  was  set  equal  to  zero  below  10'  Muv. 
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Tabic  X 


TF!ST  5  r.  OAA  PA  DfC  r\  TT  Ttm  r*  A  T,  T  T-x  r  r  _ 

A  ~  - *  ’  A  A1X-'  1  A  JVU1  i  1  r.K/v  1  1U1N 

ABSOLUTE  INTEGRAL  FLUX  RESULTS  WITH 
THE  TRUE  TEST  SPECTRUM 


N  eutron ' 
Energy 
(Mev) 

Integral  Neutron  Flux 

~  TV. .  = - 

T rue  Test 
Spectrum 
(TTS) 

SAND  II 
Spect  rum 

Ratio^ 

SAND  11  /TTS 

4  x  10"7 

_ 

1.112  x  105 

1  x  10'4 

- 

1.111  x  105 

1  x  10 

9.846  x  104 

1.04  1  x  105 

1.057 

1  x  10 

6.501  x  1 04 

8.784  x  104 

1 .3  5  2 

1  x  10"  1 

5.358  x  104 

6.678  x  104 

1.247 

1.0 

2.803  x  1C4 

4.372  x  104 

1.560 

1.5 

2.089  x  104 

2.999  x  104 

1.437 

2.0 

1.519  x  104 

1.865  x  104 

1.228 

3.0 

9.165  x  103 

9.558  x  103 

1.043 

4.0 

6.516  x  103 

7.386  x  103 

1.134 

5.0 

4.734  x  103 

6.606  x  103 

1.395 

6.0 

3  840  x  103 

6.300  x  103 

1.641 

7.0 

3.336  x  103 

5.789  x  103 

1.736 

8.0 

2.981  x  103 

4.307  x  103 

1.445 

9.0 

2.737  x  103 

2.045  x  103 

0.747 

10.0 

2.492  x  103 

8.592  x  102 

0.345 

11.0 

2.264  x  103 

7.914  x  102 

0.350 

12.0 

1.977  x  103 

7.896  x  102 

0.400 

13.0 

1.582  x  103 

7.866  x  102 

0.497 

14.0 

7.509  x  102 

6.378  x  102 

0.907 

15,0 

0 

1.070  x  102 

- 

+  ,p,0t  a11  of  the  "true"  spectrum's  energy  points  are  listed. 

'  T’huiaSSri\ined  errors  in  individual  foil  "measurements"  (listed  in 
able  III)  and  energy  regions  of  greatest  sensitivity  (listed  in 
Table  VI)  are  not  exactly  consistent  with  the  regions  of  maximum 
error  in  the  solution  spectrum.  This  is  because  the  inordinately 
large  assigned  errors  resulted  in  substantial  shifts  in  energy  re¬ 
gions  of  sensitivity,  which  in  effect  created  large  flux  depressions 
and  peaks  with  resultant  "gaps"  in  the  foil  energy  coverage.  Large 
errors  in  measured  activity  have  essentially  the  same  kind  of  ef¬ 
fect  anthmetmally  as  a  poor  initial  approximation;  the  phenomenon 
ot  shifts  m  energy  regions  of  sensitivity  caused  by  a  poor  initial 
approximation  is  discussed  in  Volume  II,  Set  turn  II-D-3 
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i  able  At 


TEST  6  -  COMPARISON  OF  THE  SAND  II  30th  ITERATION 
ABSOLUTE  INTEGRAL  FLUX  RESULTS  WITH 
THE  TRUE  TEST  SPECTRUM 


Neutron 

Energy 

(Mev) 

Integral  Neutron  Flux 

Ratio 

SAND  11/ TTS 

T rue  T est 
Spectrum 
(TTS ) 

SAND  II 
Spectrum 

4  x  10“7 

3,238  x  104 

2.914  x  104 

0.900 

1  x  10-6 

3.150  x  104 

2.771  x  104 

0.880 

-5 

4 

4 

lx  10 

2.882  x  10 

2.419  x  10 

0.839 

-4 

1  x  10 

2.547  x  104 

2.075  x  104 

0.815 

-3 

4 

4 

1  x  10 

2.127  x  10 

1.732  x  10 

0.814 

4  x  10"3 

1.814  x  104 

1.515  x  104 

0.835 

1  x  10-2 

1.600  x  104 

1.355  x  104 

0.847 

4  x  10"2 

1.351  x  104 

1.108  x  104 

0.820 

1  x  10"1 

1.192  x  104 

9.456  x  103 

0.793 

4  x  10”1 

8.324  x  103 

7.020  x  103 

0.843 

1.0 

4.054  x  103 

5.036  x  103 

1,242 

1.5 

3.307  x  103 

4.048  x  103 

1.224 

2.0 

3.117  x  103 

3.274  x  103 

1.050 

3.0 

2.204  x  103 

2.273  x  103 

1.031 

4.0 

2.013  x  103 

1.480  x  103 

0.735 

5.0 

1.466  x  103 

7.566  x  102 

0.516 

6.0 

1.180  x  103 

6.642  x  102 

0.563 

7.0 

9.184  x  102 

6.360  x  102 

0.693 

8.0 

7,044  x  102 

6.144  x  102 

0.872 

9.0 

6.331  x  102 

5.861  x  102 

0.926 

10.0 

5.756  x  102 

5.428  x  102 

0.943 

1 1.0 

5.1  10  x  102 

4.677  x  102 

0.915 

12.0 

4.412  x  102 

3.489  x  10Z 

0.791 

13,0 

3.135  x  102 

2.096  x  102 

0.669 

14.0 

7.580  x  101 

1.260  x  102 

1.662 

15.0 

0 

7.464  x  10  1 

- 

Not  all  of  the  "true"  spectrum's  energy  points  are  listed. 


for  Test  3 


27 

subsequent  SAND  II  run  for  Test  6,  the  A1  (n,p)  reaction  was  removed  and  a 
new  integral  solution  was  obtained,  which  had  deviations  more  consistent  with 
the  errors  listed  in  Table  III.  The  low  ratios  of  the  SAND  II  to  true  integral  flux 

values,  shown  in  Table  XI  between  4  and  8  Mev,  are  thus  due  mainly  to  a  large 

27  !|! 

inadvertent  error  in  the  value  of  activity  for  the  A1  (n,  p)  reaction,  beyond  the 
intentional  error  listed  in  Table  III. 

The  detector  results  given  in  Tables  V  and  VI,  and  the  integral  flux  re¬ 
sults  in  Tables  IX  and  X  for  Tests  4  and  5  respectively,  are  analogous  in  most 
respects  to  those  of  Tests  3  and  6,  and  require  no  further  discussion. 

Tests  1  through  13,  as  well  as  others  not  reported  here,  have  helped  to 
substantiate  the  expected  capability  of  the  SAND  II  code  to  yield  integral  flux  re¬ 
sults  (by  about  the  5th  to  30th  cycle  of  iteration)  that  are  accurate  to  within  10 
to  30%  at  all  energy  points  over  the  range  of  10  ^  to  18  Mev,  if  the  absolute 
errors  in  measured  activity  and  foil  reaction  cross  sections  are  within  similar 
limits  and  the  set  of  foil  reactions  has  been  properly  selected  for  the  spectrum 
being  measured. 

SECTION  V 

INTERPRETATION  OF  EXPERIMENTAL  DATA 
USING  THE  SAND  II  CODE 

A.  Foil  Activation  Measurements 
1,  Foil  Set  Selection 

For  the  tests  discussed  in  Section  V-B,  the  choice  of  foil  material  was 
dictated  primarily  by  the  half-life  of  the  product  nuclide  and  the  accuracy  with 
which  the  reaction  cross  section  was  known.  Another  important  factor  was  the 
degree  of  gamma  heating  and/or  fission  heating  which  the  foil  would  experience. 
Table  XII  is  a  compilation  of  all  the  foil  detectors  used,  along  with  pertinent 
nuclear  data. 

In  general,  all  non-fissionable  foils  were  made  ~l/2  inch  diameter  and 
then  were  packaged  together  so  that  the  axial  length  of  one  completely  assembled 

>:<ln  reference  3,  Figure  3,  the  negative  displacement  of  the  integral  flux  point 
for  the  Al27(n,p)  reaction  is  also  consistent  with  this  conclusion,  which  had  not 
been  noted  in  the  previous  study. 
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Table  XII 

FOIL  DETECTOR  NUCLEAR  AND  PHYSICAL  PROPERTIES 


Foil  Form 

Foil  Reaction 

Product 
Nuclide'  h 
Half- Life 

u238 

u238  (n,y)U239 
u238  (n.f)  F.  P. 

* 

Th232 

Th232  (n,f)  F,  P. 

s|( 

U235 

U235  (n,  f)  F.  P. 

* 

Np237 

Np237  (n,f)  F.  P. 

* 

NaCl  (30  to  187  mil  foils) 

Na23  (n,y)  Na24 

15.00  hours 

Cl35  <n,i*)  P32 

14.41  days 

Au  (1/2  and  2  mil  foils) 

Au197(n,y,Au198 

2,698  days 

0.130%  Au-Al  (20  mil  wire) 

Au197  <„,y>  Au198 

2,698  days 

Co  (40  mil  wire) 

Co59  (n, > )  Co60 

5,242  years 

n.483%  Co-Al  (30  mil  wire) 

Co59  (n,  y )  Co60 

5.242  years 

0,457%  Co-Al  (5  mil  foil) 

Co59  (n,y)  Co60 

5,242  years 

80%  Mn  -  20%  Cu  (2  mil  foil) 

Mn55  (n,  y )  Mn56 

2.58  hours 

Cu63  (n,y)  Cu64 

1 2.  84  hours 

Cu  (0,7  to  5  mil  foils) 

Cu63  (n,y)  Cu64 

1  2. 84  hours 

In  ( 1  0  mil  foil) 

t  115,  t-1 1  5m 

In  (n,  n  )  In 

4. 5 1  hours 

Fe  (5  mil  foil  and  £  20  mil  wire) 

Fe54  (n,  p)  Mn54 

313.5  days 

Fe 38  (n,  p)  Mn38 

2.58  hours 

S  ()  87  mil  foil) 

S32  (n,  p)  P32 

14,41  days 

Ni  (10  mil  foil) 

Ni58  (n,  p)  Co58 

71,7  days 

Ni58  (n,  2n)  Ni37 

36. 7  hours 

Ai  (7  and  30  mil  foils) 

Al27  (n,a)  Na24 

1  5. 00  hours 

Mg  (5  mil  foil) 

Mn  (n,  p)  Na 

15,00  hours 

K1  (70  to  150  mil  foils) 

I127(n,2n)I128 

13.1  days 

NH4I  (187  mil  foil) 

I127(n,2n)I128 

1 3. 1  days 

Zr  (5  mil  foil) 

Zr90  (n.  2n)  Zr89 

78.4  hours 

,:‘These  reactions  were  analyzed  by  fission  product  (Mo99)  and/or  fission 
track  counting  (except  for  the  U238  (n ,y)  reaction  which  was  analyzed  bv 
Np“39  counting)  by  Armani  (Ref.  18)  and  Tochilin  (Ref.  19)  respectively. 
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set  was  about  1/2  to  1  inch.  The  individual  foils  were  separated  by  1.3  mil 
aluminum  spacers  to  prevent  cross  contamination  during  handling.  The  entire 
foil  set  was  wrapped  in  aluminum  foil  and  irradiated  bare,  cadmium-covered, 
and/ or  boron-covered  as  required.  In  some  instances,  such  as  for  the  reso¬ 
nance  detectors,  individual  foils  were  prepared  and  irradiated  separately  in 
order  to  study  the  self- shielding  effect  of  the  entire  package  (or  a  specific  foil 
material).  In  this  regard  the  use  of  thick,  thin,  and  alloyed  resonance  foil  ma¬ 
terials  helped  considerably  in  the  analysis  of  data.  (Experience  gained  from 
these  tests  has  shown  that  it  is  not  practical  to  use  alloys  and/or  backing  mate¬ 
rial  in  cases  where  the  diluent  and/or  backing  materials  have  reaction  product 
nuclides  longer  lived  than  the  reaction  product  nuclide  of  interest,  such  as  indium 
with  aluminum  backing.)  The  fission  detectors  were  always  contained  in  sepa¬ 
rate  packages  that  were  placed  adjacent  to  (or  in  the  same  location  as,  if  separate 
irradiations  were  made)  the  non-fissionable  foils.  In  some  cases  involving  low 
flux  levels,  it  was  necessary  to  wrap  the  iron  and  cobalt-aluminum  wires  around 
the  outside  of  the  assembled  package  to  increase  the  total  quantity  of  foil  mate¬ 
rial  without  causing  serious  self-absorption-scattering  effects. 

An  initial  irradiation  using  most  of  the  foil  materials  uncovered  some 
problems  with  the  foil  detectors  that  might  not  have  become  evident  in  subsequent 
tests.  Examination  of  these  foils  revealed  that  cross  contamination  of  the  foil 
materials  had  occurred  during  manufacture  and  that  some  of  the  certified  assay 
values  for  the  alloys  were  inaccurate.  The  correct  values  for  these  alloys  de¬ 
termined  by  subsequent  chemical  analyses  are  given  in  Table  XII, 

»!< 

2.  Foil  Activity  Measurement  Techniques 

Comparable  capability  for  measuring  the  foil  activities  exists  at  each  of 
the  laboratory  sites  [ANL  (Ref.  18),  USNRDL  (Ref.  19),  Sandia  (Ref,  20),  BNW 
(Ref.  21),  and  Stanford  (Ref.  22)]  involved  in  this  program,  as  evidenced  by  the 
excellent  agreement  of  results  obtained  for  the  measurement  of  foils  common  to 
most  laboratories.  Each  laboratory  measured  particular  foils  or  foil  sets. 

Most  of  the  sets  were  measured  by  Atomics  International. 


*Only  the  methods  employed  at  Atomics  International  arc  discussed,  but  tln-y 
are  representative,  in  a  general  sense,  of  currently  used  techniques. 
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A  477  gamma  ionization  chamber  is  the  primary  instrument  at  AI  for  the 
calibration  of  gamma -emitting  nuclides.  Cross  calibration  of  the  chamber’s 
response  with  the  National  Bureau  of  Standards  has  been  effected  for  all  avail- 

^  2  2  (y 

able  nuclides.  The  daily  response  of  the  4lt  GIC  is  checked  against  an  NBSRa" 
standard.  Table  XIII  reports  results  of  an  interlaboratory  comparison  con¬ 
ducted  by  the  ASTM  Committee  E-10,  Sub  V,  on  Neutron  Dosimetry  (Ref.  23). 
Similar  results  have  been  reported  by  Yoshikawa  and  Zimmer  (Ref.  24). 


Table  XIII 

ASTM  INTERLABORATORY  CALIBRATION  RESULTS  FOR 
SELECTED  FOIL  REACTIONS* 


Reaction  Product 

Laboratory 

-  58 

Co 

Co60 

Mn54 

Argonne  National  Laboratory 

0.996 

0.989 

0.983 

0.979 

1.024 

1.023 

IIT  Research  Institute 

0.960J 

0.957* 

0.996 

0.971 

1.037 

0.999 

Idaho  -  NRT 

0.986 

1.000 

1.000 

0.996 

1.034 

1.004 

Atomics  International 

0.982 

0.982 

1.000 

1.017 

0.997 

1.000 

U.  S.  Naval  Radiological  Defense 
Laboratory 

1.015 

1.015 

1.015 

1.010 

1.024 

1.014 

Pratt  and  Whitney 

1.000 

0.993 

0.922 

0.917 

1.019 

0.992 

Euratom 

1.04 

(0,997)§ 

0.996 

1,05 

(1.003)§ 

0.994 

>!<The  values  given  are  the  ratio  of  the  reported  activity  by  the  participating 
laboratories  to  that  measured  by  the  reference  laboratory  for  two  separate 
foils;  Battelle  Northwest  (R.  Dahl)  served  as  the  reference  laboratory  and 
provided  the  irradiated  foils  which  had  been  measured  by  W.  Zimmer  at 
Iso-Chem  Corporation. 

fThese  values  were  obtained  assuming  a  decay  scheme  somewhat  different 
from  that  used  by  the  other  laboratories;  if  the  same  decay  scheme  had 
been  used,  these  values  would  have  been  about  2%  higher, 

KRerneasured 
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Normally  the  set  of  calibrated  gamma  spectrometer  standards  is  main¬ 
tained  at  appr uximalel y  40  nuclides'  with  many  shorter  lived  nuclides  being  pro¬ 
duced  and  calibrated  as  needed.  These  standards  serve  as  both  gamma  energy 
and  quantitative  activity  standards  for  the  spectrometer. 

At  AI,  a  477  gamma  re-entrant  high  pressure  (40  atmospheres)  ionization 
chamber,  a  ZTT  proportional  flow  counter,  and  a  400  channel  gamma  spectrom¬ 
eter  (with  associated  readout  equipment)  utilizing  both  a  3  in,  x  3  in.  Nal(Tl) 
well  crystal  and  a  3  in.  x  3  in.  Nal(Tl)  end  window  crystal,  were  involved  in 
making  the  measurements  for  these  tests. 

Since  some  of  the  foil  detectors  were  alloys  (for  reasons  discussed  in 

Section  V -A- 1 ),  the  gamma  spectrometer  was  used  to  measure  the  separate 

activities  of  the  components.  The  half-lives  were  followed  to  ascertain  that  no 

impurities  were  present  to  contribute  errors  to  the  measurement.  In  some 

foils  the  reaction  product  nuclide's  half-life  was  long  and  did  not  afford  a  means 

with  which  to  check  the  nuclide  purity.  The  foil  material  Co-Al  (see  Table  XII), 

24 

after  a  reasonable  decay  period,  was  essentially  free  of  the  Na  activity  from 
the  Al  (n,a)  reaction  and  was  chemically  free  of  other  starting  material  impuri¬ 
ties,  yet  was  still  nuclidically  impure.  Cobalt-59,  which  is  1  00%  abundant,  still 
has  three  reactions,  Co^  (n,y)Co^;  Co"^  (n,  p)  Fe^,  and  Co"^  (n,  2n)  Co^; 
the  product  nuclides'  gamma  ray  energies  from  the  last  two  are  similar  enough 
to  interfere  with  the  1.17  Mev  and  1,33  Mev  photopeaks  commonly  used  for  the 
measurement  of  Co^.  To  avoid  errors  caused  by  either  over-  or  under¬ 
corrections  for  these  extraneous  activities,  the  2.50  Mev  coincidence  sum  peak, 
which  is  well  removed  from  interference  from  these  nuclides,  was  utilized  for 
some  of  the  measurements. 

Selection  of  Zn  as  a  foil  detector  posed  a  similar  problem  due  to  the 

y-energy  and  half-life  similarities  of  the  reaction  product  nuclides  from  the 

Zn^  (n,  p)  Cu^  and  (n,y)  Zn^m  reactions.  Zinc  69m  emits  a  gamma  ray 

64 

with  an  energy  of  0,44  Mev  and  has  a  half-life  of  13.9  hours,  whereas  Cu  emits 

a  positron  and  has  a  half-life  of  12,8  hours.  (The  1,34  Mev  gamma  ray  from 
64 

Cu  is  not  useful  because  of  the  relatively  low  number  of  gammas  per  disinte- 

64 

gration.)  When  the  positron  from  Cu  is  annihilated  two  gamma  photons  of 
0,51  Mevare  emitted.  By  using  a  NaI(T  1 ) well  crystal,  the  sum  peak  at  1.02  Mev 
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due  ro  me  gamma  photon  pair  from  the  positron  annihilation  can  be  used  in  the 
64 

determination  of  Cu  activity,  thus  avoiding  the  interference  at  0.5  1  Kiev  caused 
by  the  0.44  Muv  gamma  ray  from  Zn 

This  technique  involving  the  measurement  of  positron  annihilation  radiation 
must  be  used  with  care  to  assure  that  no  other  positron-emitting  nuclide  is  pres¬ 
ent  to  interfere.  Such  interference  may  be  difficult  to  resolve,  especially  if  the 

•  fy  s 

half -lives  of  the  two  positron-emitters  are  similar.  In  zinc,  the  Zn  (n,y)Zn 

produces  a  positron-emitting  nuclide,  although  in  this  case  it  may  be  corrected 

64 

for  easily  because  of  the  very  different  half-lives  (12.9  hours  for  Cu  and  243 
days  for  Zn^"*). 


In  certain  foils,  it  may  be  impossible  to  distinguish  an  err  or -contributing 

impurity  because  a  reaction  product  from  the  impurity  is  the  same  as  the  product 

56  56 

of  the  pertinent  reaction.  For  example,  the  Fe  (n,  p)Mn  activity  measured  in 

55  56 

an  iron  foil  may  have  interference  from  the  Mn  "(n,y)Mn  reaction  caused  by 

the  presence  of  manganese  as  an  impurity;  in  this  case,  manganese  quantities  of 

56 

only  several  hundredths  of  one  percent  may  yield  Mn  of  the  same  order  of  mag¬ 
nitude  as  that  produced  by  the  reaction  of  interest.  In  one  of  the  experiments 
discussed  in  Section  V-B-2  (ECEL  test),  this  problem  was  made  evident  by  a 

comparison  of  the  results  of  bare  foil  and  boron  10-covered  foil  irradiations;  the 
56 

measured  Mn"  activity  was  a  factor  of  two  higher  for  the  bare  iron  foil  than 

55  56 

for  the  one  in  a  boron  10  cover.  (The  Mn  (n,  y)Mn  reaction  cross  section  is 
most  significant  in  the  resonance  region  and  below;  flux  in  this  region  was  es¬ 
sentially  removed  by  the  boron  10  cover,  ) 

In  some  tests  the  high  cross  section  of  In*  *  "*(n,  y)In*  *°m  necessitated  a 
24-hour  cooling  period  before  attempting  the  In**"’^,  n /)In‘  determination. 
Indium  1.14  produced  by  the  In*  *  ^(n,  y)Ts*  *^m  reaction,  because  of  its  relatively 
long  half-life  and  the  ~-5  half-life  delay  imposed  by  In**km}  was  approximately 
an  order  of  magnitude  greater  in  activity  than  the  In*  This  three-way  de¬ 

cay  permitted  approximately  an  8-hour  period  (starting  24  hours  after  time  of 
removal)  in  which  to  determine,  the  In**"’111  activity  and  to  verify  the  half-life. 

In  spectra  that  were  harder,  a  delay  of  10  hours  eliminated  the  In**^m  inter - 

„  ,,  . .  ,  i  1 5 .  Mt  115rn  .  t  .  .  ,  ,  T  114m 

ference,  allowing  the  In  (n,  n  )ln  measurement  to  be  made  before  In 

interfered  appreciably. 
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Pure  beta -emitting  nuclides  were  determined  by  measuring  either  n  sub¬ 
limed  or  evaporated  aliquot  (on  a  platinum  disk)  in  a  2ff  proportional  flow  boom 
ter.  Other  means  of  measurement  exist,  but  this  technique  is  quite  adequate 
for  beta  energies  greater  than  ~  0.3  Mev.  As  in  the  previous  cases  discussed, 

nuclide  purity  is  a  concern  and  now  a  much  more  serious  problem  than  it  had 

35  32 

been  in  y-spectrometric  measurements.  For  example,  the  Cl  (n,  a)P  reac- 

35  35  32  32 

tion  has  interference  from  Cl  (n,  p)S  ,  and  the  S  (n,  p)P  reaction  has  in- 
33  33 

terference  from  S  (n,  p)P  ;  all  four  product  nuclides  are  pure  beta  emitters. 
32 

The  P  activity  was  determined  in  both  cases  using  an  aluminum  absorber  to 

35 

remove  the  contribution  of  the  weaker  betas  (0.168  Mev  from  S  and  0.25  Mev 
33 

from  P  )  to  the  count  rate,  and  subsequently  correcting  the  count  rate  for  the 
32 

attenuation  of  the  P  betas  (1.71  Mev)  by  the  absorber. 

This  discussion  does  not  exhaust  the  possible  difficulties  that  might  be 
encountered  in  foil  activity  measurement.  It  is  clear  that,  when  reactions  are 
used  with  which  there  is  little  or  no  previous  experience,  or  when  neutron  en¬ 
vironments  are  encountered  which  might  produce  unforeseen  effects,  the  mate¬ 
rials  used  and  the  reactions  measured  must  be.  carefully  examined  to  determine 
whether  any  problems  such  as  those  discussed  here  may  be  in  effect,  before  any 
conclusions  are  drawn  regarding  the  accuracy  of  the  actual  activity  measure¬ 
ments  or  cross  section  evaluations. 

It  is  interesting  from  a  practical  point  of  view  to  observe  that  the  type 
of  agreement  achieved  in  the  ASTM  interlaboratory  comparison,  Table  XIII, 
was  never  achieved  in  the  current  tests.  There  are  reasons  for  this,  however; 
the  ASTM  study  only  involved  the  re-counting  of  foils  that  had  been  previously 
counted  by  the  reference  laboratory  (BNW).  For  the  tests  discussed  in  this  re¬ 
port,  Section  V-B,  the  results  are  based  on  the  counting  of  different  foils  by 
different  laboratories;  i.  e.  ,  each  laboratory  provided  its  own  sealed  foil  set 
that  was  irradiated  and  returned  for  counting.  Since  multiple  irradiations  were 
also  performed  in  some  cases,  the  assigned  values  of  error  reflected  errors 
due  to  changes  in  flux  level  and  exposure  times.  Consequently  the  values  of  er¬ 
rors  (one  standard  deviation)  listed  in  Tables  XIV,  XVI,  XVIII,  XX,  XXI,  XXII, 
XXIII,  XXIV,  and  XXV  are  realistic,  in  the  sense  that  they  indicate  what  cur¬ 
rently  can  be  achieved  (with  a  reasonable  effort)  in  cases  where  the  multiple 
foil  technique  has  been  utilized. 
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B.  SAND  II  Foil  Activation  Data  Reduction 


1.  Summary  Discussion 

The  SAND  II  code  has  been  used  to  determine  experimentally  the.  neutron 
flux  spectra  in  different  types  of  neutron  environment.  In  each  environment,  a 
set  of  foil  detectors,  utilizing  both  low  and  high  energy  neutron  reactions,  was 
irradiated.  Measured  foil  activities  for  each  foil  set  (extrapolated  to  saturation, 
except  in  one  pulsed  case,  in  which  time-of-removal  values  were  required)  were 
used  as  input  to  the  code  to  obtain  an  appropriate  solution  flux  spectrum  for 
each  environment. 

One  of  tlie  goals  of  the  present  work  was  lo  compare  the  results  of  the 
foil  activation  method  with  spectrometer  measurements,  muitigroup  diffusion 
and  transport  calculations,  and  Monte  Carlo  calculations.  Emphasis  was  placed 
on  attempting  to  correlate  the  effect  of  scattering  and  absorption  (in  materials 
such  as  carbon,  oxygen,  beryllium,  iron,  aluminum,  ch  rnium,  and  other  ele¬ 
ments  present  in  a  given  environment)  with  structure  that  appeared  in  the  SAND  II 
iterative  solution  spectra.  This  required  that  a  distinction  be  made  between 
possible  "true"  structure  and  that  resulting  from  foil  activity,  foil  self-absorption, 
and  foil  cross  section  measurement  errors. 

As  a  convenience  in  the  analysis  of  the  experimental  data,  the  thirtieth 
iterative  results  are  used  as  the  solution  for  the  first  five  cases  studied.  For 
the  sixth  case,  SAND  II  results  are  reported  that  were  achieved  using  solution 
and  foil  discard  criteria  of  10%  and  1.96  (95%  confidence  level),  respectively. 

For  the  experimental  cases  studies,  the  foil  irradiations  were  performed 
using  tlie  following  facilities: 

a.  WMGR,  at  the  center  of  a  homogeneous  Water  Moderated  and  Graphite 
Reflected  water  cooled  reactor; 

b.  GMRA,  at  the  center  of  a  large  graphite  central  region  in  a  Graphite 
Moderated  and  Reflected  Air  cooled  reactor; 

c.  GMRW,  in  a  test  location  in  a  Graphite  Moderated  and  Reflected 
Water  cooled  reactor; 

d.  SMFS,  in  a  collimated  beam  from  a  Slightly  Moderated  Fission 
Source ; 
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e.  BRGM,  in  the  leakage  spectrum  from  a  Beryllium  oxide  Reflected 
Graphite  Moderated  gas  cooled  reactor; 

f.  ECEL,  in  the  central  core  region  of  the  Epithermal  Critical  Experi¬ 
ment  Laboratory  fast  critical  assembly,  Core  14-13. 

In  the  SMFS  beam,  two  tests  were  conducted  in  somewhat  different  loca¬ 
tions  (expected  to  yield  qualitatively  similar  results,  except  for  differences  in 
side-  and  back-scattering  effects).  In  both  the  GMRW  and  BRGM  tests,  there 
were  significant  amounts  of  foreign  materials;  in  the  GMRW,  the  foils  used 
were  surrounded  by  approximately  1/4  inch  each  of  stainless  steel,  aluminum 
and  water,  and  in  the  BRGM  the  foils  were  separated  from  the  source  by  approx¬ 
imately  1/2  inch  of  Hastelloy  C.  The  core  region  of  the  ECEL  was  made  up  of 
repeating  2"  x  2"  (by  1/8  to  1/4"  thick)  squares  of  carbon,  aluminum,  stain¬ 
less  steel,  uranium  238,  and  uranium  235  (Ref.  25).  (Because  of  this  arrange¬ 
ment  of  materials,  the  ECEL  core  region  could  be  considered  to  be  homogeneous 
for  reactor  physics  calculations.  ) 

It  is  clear  from  the  study  and  analysis  of  these  six  test  cases  that  most 
materials  present  in  a  given  neutron  environment  have  certain  identifiable  ef¬ 
fects  on  the  flux  —  that  the  total  cross  sections  of  all  the  materials  in  an  envi¬ 
ronment  are,  in  a  rough  qualitative  way,  mirrored  in  the  neutron  spectrum. 

Some  examples  of  this  may  be  seen  in  Fig.  11.  All  the  structural  features  of 
the  calculated  spectra  shown,  as  well  as  most  of  those  of  the  measured  spectra, 
are  quite  readily  identifiable  with  structure  in  the  total  cross  sections  of  the 
materials  known  to  be  present.  The  spectra  shown  in  Fig.  11  will  be  used  for 
comparison  in  discussing  the  test  results  presented  in  Section  V-B-2. 

2.  Iterative  Results 

The  SAND  11  foil  activation  results,  integral  flux  comparison,  and  dif¬ 
ferential  flux  comparison  for  the  WMGR,  GMRA,  and  GMRW  tests  are  presented 
in;  Tables  XIV,  XVI,  and  XVIII;  Tables  XV,  XVII,  XIX;  and  Figs.  12,  13,  and 
14,  respectively.  In  each  figure,  the  dashed  curve  is  the  theoretical  calculation 
(used  as  the  input  approximation)  and  the  solid  curve  is  the  30th  iteration  solu¬ 
tion.  In  the  lower  part  of  each  curve,  pertinent  representative  spectra  shown 
in  Fig.  11  are  reproduced  (as  dotted  circle  curves)  to  help  establish  the  credi¬ 
bility  of  some  of  the  differential  structure  generated  in  the  iterative  solution. 
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HETEROGENEOUS  LIGHT  WATER  MODERATED  AND 
COOLED  THERMAL  REACTOR  i  PROTON  RECOIL  EMULSION 
(Ref.  ?Si]  • 


DEEP  PENETRATION  IN  H  Q  I  CALCULATION  (Ref.  27)1 


TRIGA  MARK  III  LEAKAGE  SPECTRUM, 

2  in,  H20  REFLECTOR  [SPECTROMETRY  (Ref.  15)] 

SWIMMING  POOL  REACTOR.  CLOSE  TO  CORE. 

[SPECTROMETRY  (Ref.  28)] 

DEEP  PENETRATION  IN  HYDROGEN  [CALCULATION  (Ref.  27)] 


DEEP  PENETRATION  IN  <CH)x  [CALCULATION  (Ref,  27)1 
SCATTERED  BY  CARBON  AT  70°  [SPECTROMETRY  (Ref,  29)]  • 

DEEP  PENETRATION  IN  CARBON  [CALCULATION  (Ref.  27)]  • 

DEEP  PENETRATION  IN  BERYLLIUM  [CALCULATION  (Ref.  27)]  • 

SCATTERED  BY  ALUMINUM  AT  70°  [SPECTROMETRY  (Ref,  29)]  • 
SCATTERED  BY  STEEl  AT  i  11' '  SPECTROMETRY  (Ref.  29)1  • 


Figure  11,  Representative  Neutron  Flux  Spectra  Tor  Different  Envi  mnnu’iil  s , 
Showing  Depression  Effects  of  Material  Cross  Sections 
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■  ARBJTRARY  SCA'lEi 


Table  XIV 
WMGR  TEST 


SAND  II  THIRTIETH  ITERATION  RATIOS  OF  MEASURED  TO 
CALCULATED  ACTIVITIES 


Foil(a* 

Reaction 

Saturated 

Measured 

Activity 

(Dps/Nucleus) 

Saturated 
Calculated 
Activity 
(Dps/Nurleus ) 

Ratio 

Meas,  /Calc. 
Activities 

Na23 

(n,  y)  Na24 

4.21 

X 

10- 26 

± 

5% 

4.233 

X 

10_ 

26 

0.995 

Mn55 

(n,  y)  Mn^ 

2.02 

X 

10"  24 

± 

7% 

2.031 

X 

.10" 

24 

0.995 

Cu63 

(n,y)  Cu64 

7.12 

X 

IQ*25 

± 

12% 

7.158 

X 

10" 

25 

0.995 

Co59 

(n,y)  Co^° 

1.01 

X 

10”23 

± 

3% 

1.015 

X 

10 

23 

Au19' 

7(n,y)Au198 

2.20 

X 

io"22 

± 

10% 

2.212 

X 

10 

22 

0.995 

Au11! 

5  ,  t  \  *11 5m 

(n, n  )  In 

2.65 

X 

10-25 

± 

10% 

2.763 

X 

10 

25 

0.959 

Fe54 

54 

(n,  p)  Mn 

2.04 

X 

io-25 

± 

15% 

1.687 

X 

10 

25 

1.21 

Ni*8 

(n,  p)  Co58 

2.14 

X 

10-25 

± 

3% 

2.012 

X 

10 

25 

1.06 

S32  (, 

*,P)P32 

7.50 

X 

IO"26 

± 

5% 

1.011 

X 

10 

25 

0.742 

Cl35 

<nf  a)  P32 

2.92 

X 

10“2fc 

± 

5% 

2.747 

X 

10 

26 

1.06 

Al27 

(n,  Of)  Na24 

1.36 

X 

IO'27 

± 

5% 

1.254 

X 

10 

27 

1.09 

24 

Mg 

(n,  p)  Na24 

2.60 

X 

io-27 

± 

5% 

2.834 

X 

10 

27 

0.918 

^127 

(n,  2n)  I126 

1.94 

X 

io-27 

± 

10% 

1.964 

X 

10 

27 

Standard  Percent  Deviation  of  the  Ratios  =  10.6 


(a)  The  foils  were  Irradiated  inside  a  ~  40  mil  cadmium-covered 
capsule. 

(b)  All  measurements  made  at  Atomics  International.  Experimental 
errors  listed  are  one  standard  deviation. 
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Table  XV 
WMGR  TEST 

SAND  II  THIRTIETH  ITERATION  NORMALIZED  INTEGRAL 
NEUTRON  FLUX  COMPARISON 


Energy  (Mev) 

Theoretical  (Ref.  30) 

Experimental 
(30th  Iteration) 

Ratio 

Exp.  /Theo, 

4.14  x  10'7 

1.000 

1.000 

mm 

1.13  x  10'6 

0.968 

0.961 

3.06  x  10"6 

0.936 

0.920 

wfim 

8.32  x  10‘6 

0.903 

0.881 

0.98 

2.26  x  10“5 

0.869 

0.840 

0.97 

6.14  x  10'5 

0.835 

0.801 

0.96 

4.54  x  10"4 

0.763 

0.725 

0.95 

3.36  x  10“3 

0.690 

0.644 

0.93 

1.7  x  l(f2 

0.626 

0.572 

0.91 

1.0  x  10“  1 

0.540 

0.475 

0.88 

4.0  x  10“  1 

0.420 

0.348 

0.83 

9.0  x  10“  1 

0.302 

0.282 

0.93 

1.  4 

0.223 

0.250 

1.12 

3.0 

7.54  x  10“2 

9.33  x  10“2 

1.24 

3.85  x  10“2 

6,37  x  10~2 

1.65 

5.0 

1.91  x  10“2 

3.84  x  10“2 

2.01 

6.0 

9.23  x  10“3 

1.32  x  10“2 

1.43 

7.0 

4.38  x  10“ 3 

5.96  x  10“3 

1.36 

8.0 

2.05  x  10“3 

2.85  x  10'3 

1.39 

9.0 

9.47  x  10“4 

1.66  x  10*3 

1.75 

-4 

,  -4 

10.0 

4.32  x  10  4 

8.60  x  10 

1.99 

11.0 

1.95  x  10“4 

3.78  x  10'4 

1.94 

12.0 

8.72  x  10“5 

1.83  x  10“4 

2.09 

13.0 

3.84  x  10“  5 

8.11  x  H)“6 

2. 1 1 

14.0 

1.66  x  10“5 

3.57  x  10“5 

2. 15 

15.0 

6.92  x  10“6 

1.51  x  lO'3 

2.  18 

16.0 

2.66  x  10“6 

5.82  x  10“6 

2.  19 

17.0 

8.04  x  10'7 

1.76  x  U)'6 

2.  18 

2- 


Table  XVI 
UMKA  TEST 

SAND  II  THIRTIETH  ITERATION  RATIOS  OF  MEASURED  TO 
CALCULATED  ACTIVITIES 


Foil(a) 

Reaction 

Saturated 
Measured 
Activity 
(Dps  /Nucleus ) 

Saturated 
Calculated 
Activity 
(Dps  /Nucleus ) 

Ratio 

Meas.  /Calc. 
Activities 

m  23 
Na 

(n,  y) 

24 

Na 

3.26 

X 

1C’27 

± 

5% 

3.276 

X 

10-27 

0.995 

Mn55 

(n,y) 

MnS6 

1.63 

X 

IQ'25 

± 

7% 

1.638 

X 

10‘25 

0.995 

Cu63 

(n,y) 

Cu64 

6.00 

X 

io’26 

± 

10% 

6.029 

X 

io-26 

0.995 

Co59 

(n.y) 

Co60 

1.15 

X 

10-24 

± 

20% 

1.156 

X 

jo"24 

0.995 

a  197,  , 

Au  (n,y) 

Au198 

1.90 

X 

io*23 

± 

7% 

1.909 

X 

10‘23 

0.995 

In115 

(n,  n  ' 

)InU5m 

1.75 

X 

10-27 

± 

10% 

1.801 

X 

io-27 

0.972 

Ni5S 

(n,  p) 

_  58 

Co 

9.62 

X 

10-28 

± 

3% 

8.004 

X 

io-28 

1.20 

s32 

(n,  p) 

p32 

2.88 

X 

10-28 

5% 

3.737 

X 

10*28 

0.771 

Cf35 

(n,  a) 

p32 

1.09 

X 

10" 28 

± 

5% 

9.985 

X 

10‘29 

1.09 

A!27 

(n,  a) 

Na24 

4.69 

X 

10-30 

± 

10% 

4.677 

X 

io'30 

1.00 

Mg24 

(n,  p) 

Na24 

1.13 

X 

10-29 

± 

10% 

1.147 

X 

io-29 

0.985 

Standard  Percent  Deviation  of  the  Ratios  =  10,2 


(a)  The  foils  were  irradiated  inside  a  60  mil  cadmium -cove red 
rotating  capsule  that  provided  identical  exposures  for  all  foils, 
which  were  located  in  a  plane  at  the  same  radial  distance. 

(b)  All  measurements  made  at  Atomics  International.  Experi¬ 
mental  errors  listed  are  one  standard  deviation. 
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Table  XVII 
GMRA  TEST 


SAND  II  THIRTIETH  ITERATION  NORMALIZED  INTEGRAL 
NEUTRON  FLUX  COMPARISON 


I 

\:v 

Energy  (Mev) 

Theoretical  (Ref.  31) 

Experimental  (b) 
(30th  Iteration) 

Ratio 

Exp.  /'Theo. 

f 

f 

4. 14 

x  10'7 

1.000 

1.000 

1.00 

Y 

I 

1. 13 

x  10'6 

0,027 

0.968 

1.04 

3.06 

x  10'6 

0.853 

0.905 

1.06 

8.32 

x  10'6 

0.780 

0.831 

1.07 

v- 

2.26 

x  10‘5 

0.706 

0.758 

1.07 

}V 

6. 14 

x  10‘5  ' 

0.635 

0.682 

1.07 

l 

4. 54 

x  10‘4 

0.494 

0.532 

1.08 

A 

3.36 

x  10“3 

0.362 

0.390 

1.08 

l 

* 

1.7  x 

10-2 

0.265 

0.273 

1.03 

?•  * 

1.0"x 

10“ 1 

0.171 

0.156 

0.91 

r 

4.0  x 

10"  1 

0.107 

7.00  x 

ID'2 

0.65 

•)  * 

9.0  x 

10“ 1 

5.12  x 

10“2 

4.27  x 

10“2 

0.83 

1.4 

3.17  x 

10‘2 

2.78  x 

10“  2 

0.88 

Ij 

; 

3.0 

3.96  x 

10“3 

6.25  x 

10“3 

1.58 

!. 

4.0 

2,02  x 

10“3 

4.37  x 

10“  3 

2.  16 

s 

5.0 

1.00  x 

10“  3 

2.62  x 

10“3 

2.62 

i 

l 

6.0 

4.85  x 

10“4 

1.00  x 

10“3 

2.06 

J- 

} 

7.0 

2.30  x 

10“4 

4.57  x 

io“4 

1.99 

\ 

8.0 

1.08  x 

10"4 

1.96  x 

10“4 

1.81 

( 

9.0 

4.97  x 

10“5 

8.55  x 

10“5 

1.72 

10.  0 

2.27  x 

-5 

10  3 

3.67  x 

10“  5 

1.62 

11.0 

1.03  x 

10-6 

1.59  x 

10“  5 

1.54 

12.0 

4.58  x 

10"6 

6.88  x 

10“6 

1.50 

1 

13.0 

1.79  x 

10“6 

2.93  x 

lo'6 

1.64 

;  • 

14.0 

8.72  x 

10“7 

1,24  x 

10"6 

1.42 

> 

15.0 

3,64  x 

10'7 

5.05  x 

10“7 

1.39 

- 

16.0 

1.40  x 

10'7 

1.89  x 

10“7 

1.35 

17.0 

4.22  x 

10“8 

5,56  x 

10"8 

1.32 
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Figure  13,  SAN  U  II  Differential  Flux  Results  for  the  Cl  MR  A  Test 


Table  XVIII 
GMRW  TEST 


SAND  II  THIRTIETH  ITERATION  RATIOS  OF  MEASURED  TO 
CALCULATED  ACTIVITIES 


Reaction 


CoJ/ 

,,238 


Saturated 
Measured 
Activity 
(Dps /Nucleus) 

4.40 

X 

10-23 

± 

6% 

2.17 

X 

10-24 

± 

12% 

9.42 

X 

IQ"24 

± 

10% 

8.84 

X 

ID'24 

± 

5% 

4.12 

X 

io"*25 

5% 

9.03 

X 

10-27 

5% 

1.36 

X 

10'25 

5% 

1.37 

X 

10_  26 

± 

5% 

2.32 

X 

10"2* 

± 

5% 

6.11 

X 

10'27 

± 

5% 

3.42 

X 

10-27 

± 

9% 

5.98 

X 

10“27 

± 

4% 

7.22 

X 

io"2* 

± 

5% 

2.91 

X 

10-29 

5% 

6.29 

X 

10“  29 

± 

5% 

6.67 

X 

10 -z9 

i: 

5% 

Saturated 
Calculated 
Activity 
(Dps/Nucleus ) 


Ratio 

Meas.  /Calc. 
Activities 


(n ,7)  Au  7 
(n,  y)  Co60 
(n,  y)  Np23 
(n,f)  F.P. ' 
(n,  y)  Mn^ 
(n,  y)  Na2' 
(n,  y)  Cu^ 


(n,  f )  F.P. 
(n,f)  F.P. 
(n,  p)  Mn5 
(n,  p)  Co5{ 


(n,  a)  Na 
(n,  p)  Na 


Standard  Percent  Deviation  of  the  Ratios  =  7.98 

(a)  The  foils  were  irradiated  inside  an  ~  40  mil  cadmium-covered 
capsule. 

(b)  Measurements  made  at  Atomics  International,  Argonne  National 
Laboratory  (Ref.  18),  Battelle  Northwest  (Ref.  21),  and  Sandia 
Corporation  (Ref.  20).  Experimental  errors  listed  are  one 
standard  deviation. 

(c)  This  reaction  was  not  used  to  obtain  the  iterative  solution  since 
a  correction  for  self-absorption  has  yet  to  be  determined 
experimentally. 


4.413 

X 

10 

2.177 

X 

10' 

9.449 

X 

10' 

4.133 

X 

10“ 

9.057 

X 

10' 

1.364 

X 

10“ 

1.411 

X 

10' 

2.372 

X 

10' 

5.441 

X 

10" 

4.368 

X 

10“ 

5.868 

X 

10" 

6.242 

X 

10“ 

2.871 

X 

10" 

6.426 

X 

10" 

6.703 

X 

10" 

Table  XIX 
GMRW  TEST 


SAND  II  THIRTIETH  ITERATION  NORMALIZED  INTEGRAL 
NEUTRON  FLUX  COMPARISON 


Energy  (Mev) 

Theoretical  (Ref.  32) 

Experimental 
(30th  Iteration) 

Exp.  /Theo. 

1.30  x  10'4 

0.779 

0.704 

0.90 

4.54  x  10-4 

0.721 

0.637 

0.88 

5.53  x  10"3 

0.605 

0.464 

0.77 

1,50  x  10"2 

0.563 

0.406 

0.72 

4.09  x  10'2 

0.518 

0.345 

0.67 

1.10  x  10-1 

0.463 

0.278 

0.60 

2.35  x  lO"1 

0.405 

0.210 

0.52 

4.98  x  10"  1 

0.328 

0.141 

0.43 

1.05 

0.221 

0.104 

0.47 

1.35 

0.184 

8.28  x  10'2 

0.45 

1.74 

0.144 

4.79  x  10*2 

0.33 

2.23 

0.105 

2.86  x  10-2 

0.27 

2.87 

6.72  x  10'2 

1.64  x  10-2 

0.24 

3.68 

4.08  x  10'2 

1.09  x  10-2 

0.27 

4.72 

2.14  x  10'2 

5.03  x  10-3 

0.24 

6.07 

7.83  x  10‘3 

1.69  x  10"3 

0.22 

8.0 

1.78  x  10'3 

4.07  x  10-4 

0.23 

9.0 

8.20  x  10'4 

2.59  x  10-4 

0.32 

10.0 

3.74  x.  10‘4 

1.72  x  10‘4 

0.46 

11.0 

1.69  x  10‘4 

0.54  x  10'5 

0.56 

12.0 

7.54  x  10'5 

4.74  x  10'5 

0.63 

13.0 

3.32  x  10'5 

2.29  x  10'5 

0.69 

14.0 

1.44  x  10"5 

1.08  x  10-5 

0.75 

15.0 

6.00  x  10"6 

4.88  x  10"6 

0.81 

16.0 

2.30  x  10'6 

2.02  x  10‘6 

0.88 

17.0 

6.96  x  10‘7 

6.53  x  10'7 

0.94 
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It  is  not  expected  that  any  of  the  spectra  from  Fig.  11  would  agree  in  struc¬ 
ture  exactly  with  the  solution  spectra  for  these  test  cases,  but  the  interactions 
among  the  effects  of  the  presence  of  water,  carbon  and  other  absorbing  /scattering 
materials  should  produce  structure  which  has  some  recognizable  characteristics. 

Some  of  the  fine  structure  detail  (superimposed  on  the  solution  envelope)  in 
Figs.  12  ,  13,  and  14  can  be  associated  with  the  properties  of  certain  foil  de¬ 
tectors.  In  Figs.  12  and  14,  for  example,  the  peak  and  dip,  respectively,  at 

~  1.7  Mev  are  known  to  be  due,  at  least  in  part,  to  errors  in  the  evaluated  cross 
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section  for  the  Fe  (n,  p)Mn  reaction.  This  type  of  effect  is  discussed  in  more 
detail  in  the  presentation  of  the  ECEL  test  results. 

The  integral  flux  comparisons  in  Tables  XV,  XVII,  and  XIX  show  rather 
substantial  disagreement  between  the  theoretical  and  SAND  II  results  over  cer¬ 
tain  energy  regions.  (In  all  three  cases,  the  results  have  been  arbitrarily  nor¬ 
malized  at  the  low  end  of  the  energy  region.  )  The  theoretical  results  would  not 
be  expected  to  be  particularly  accurate,  however,  in  any  of  these  cases,  since 
all  three  facilities  have  nuclear  and/or  structural  features  (heterogeneity)  that 
require  certain  simplifying  assumptions  to  make  the  physical  model  mathemat¬ 
ically  tractable. 

The  spectral  results  for  the  SI/IFS  and  BRGM  tests  (limited  to  the  energy 
region  between  approximately  0.1  and  10  Mev)  are  shown  by  the  solid  curves  in 
Figs.  15  and  16,  respectively.  In  the  lower  part  of  Fig.  16,  spectra  are  repro¬ 
duced  from  Fig.  11  that  show  the  type  of  flux  depressions  caused  by  carbon, 
oxygen,  and  beryllium. 

For  the  SMFS  test,  previous  time-of-flight  results  for  a  similar  neutron 
source  were  available  and  were  used  as  the  input  approximation,  represented  by 
the  dashed  curve  in  Fig.  15.  The  results  obtained  for  two  slightly  different  lo¬ 
cations  in  the  beam  with  different  side-  and  back-scattering  effects  are  shown 
and  are  identified  as  Tests  A  and  B.  The  foil  sets  and  covers  used  for  Tests  A 

and  B  were  different;  in  particular,  the  fission  foils  in  Test  A  had  a  boron  10 
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cover,  and  the  Fe  (n,p)Mn‘  reaction  was  omitted.  This  was  a  fortunate 


*The  solid  and  dashed  curves  in  each  figure  are  pertinent  to  this  discussion;  the 
upper  curves  will  be  discussed  in  Section  V-B-4. 
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Figure  15,  SAND  II  Differential  Flux  Results  lor  the  SMFS  Test 


NEUTRON  DIFFERENTIAL  FLUX  a  on-sec  Siev.f 


INPUT  SPECTRUM  7  (VOL.  IV). 


SAND  II  SOLUTION. 


REFERENCE  29,  FIGURE  11, 
SCATTERED  BY  CARBON  AT  70°. 


REFERENCE  27,  FIGURE  11,  DEEP 
PENETRATION  OF  BERYLLIUM. 


REFERENCE  26,  FIGURE  11. 
HETEROGENEOUS  LIGHTWATER 
MODERATED  AND  COOLED  THERMAL 
REACTOR. 
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Figure  16,  SAND  II  Differential  Flux  Re  suit  h  for  the  BRGM  Test 


omission;  the  most  important  difference  between  the  iterative  solutions  lor  these 

two  tests  appears  at  ~  1.7  Mev  in  the  form  of  a  spectral  peak  fur  'J  est  B,  which 

helped  to  identify  suspected  errors  in  the  evaluated  cross  sec  tion  for  the 
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Fe  (n,  p)Mn  reaction;  re-evaluation  of  the  cross  section  data  for  this  reac¬ 
tion  is  discuosed  in  Section  V-B-4. 

The  foil  reactions,  measured  and  SAND  II  calculated  values  of  activity,  and 
ratios  of  these  activities  for  Teats  A  and  B  are  given  in  Tables  XX  and  XXI,  r«- 
spactively.  Since  a  pulsed  source  was  used  for  this  test,  it  also  provides  an  ex¬ 
ample  of  the  time-integrated  mode  of  operation  for  SAND  II  using  t)ic  time-of- 
removal  (end  of  pulse)  values  of  measured  activity  as  input.  (Sec  Section  U-B). 

In  addition  to  the  tlme-of-flight  data,  comparisons  were  made  for  the  SMFS 
teat  with  the  results  of  a  Monte  Carlo  calculation  performed  by  BiggerB  (Ref.  14). 
Neither  the  time-of -flight  nor  the  Monte  Carlo  calculations  included  any  correc¬ 
tions  for  the  scattering  effect  of  adjacent  materials  In  the  SMFS  beam,  and  ex¬ 
act  agreement  between  these  and  the  Iterative  results  should  not,  therefore,  be 
expected.  The  agreement  among  these  various  integral  flux  results  (on  an  abso¬ 
lute  basis)  was  found  to  be  within  approximately  30%  over  the  energy  region 
from  10“^  to  15  Mev.  In  another  absolute  comparison  of  multiple  foil  Iterative 
results  with  Monte  Carlo  results  for  a  different  beam  type  measurement,  agree¬ 
ment  of  approximately  20%  was  found  (Ref.  33).  These  and  other  results  (in  par¬ 
ticular,  the  ECEL  test)  have  helped  to  Increase  confidence  in  the  reliability  of 
the  iterative  method  and  the  evaluated  cross  section  data  currently  In  use. 

For  the  BRGM  test,  Fig.  16,  which  involved  measurements  of  a  leakage 
spectrum,  no  reliable  information  was  available  on  the  form  of  the  unknown 
spectrum.  It  was  known,  however,  that  it  would  be  a  fission  spectrum  degraded 
by  the  graphite  moderator  and  beryllium  oxide  reflector  of  the  reactor.  A 
1/E  +  fission  form  (reference  spectrum  7,  Volume  IV)  was  therefore  specified 
as  the  input  approximation;  this  choice  was  supported  by  a  preliminary  com¬ 
parison  of  the  BRGM  measured  and  spectrum  7  calculated  values  of  saturated 
activities.  The  input  approximation  and  the  SAND  II  solution  spectrum  arc* 
shown  In  Fig.  16  by  the  dashed  and  solid  curves,  respectively;  the  foil  activity 
results  are  given  in  Table  XXII. 

The  slight  dip  in  the  solution  spectrum  at  ~  1.7  Mev  can  be  associated 
with  a  low  value  of  measured  activity  and/or  possible  cross  section  errors  for  the 
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8MFS  TEST  A 

SAND  II  THIRTIETH  ITERATION  RATIOS  OF  MEASURED  TO 
CALCULATED  ACTIVITIES 


Foil  (a) 
Reaction'  ' 

Measured^ 

Activity 

(Dps/Nucleus) 

Calculated 
Activity 
(Dpa/Nucleus ) 

Ratio 

Meas.  /Calc, 
Activities 

Na23  (n, y)  Na24 

2.10  x  10‘32  ±  15% 

2.099  x  10'32 

1.000 

AuI97(n,y)  Au198 

7.96  x  10'30  ±  10% 

7.961  x  10'30 

1.000 

Pu239(n,f)  F.P. 

3.03  x  10'24  ±  7% 

3.042  x  10'  24 

0.996 

U238  (n,  f)  F.P. 

4.66  x  10“25  ±  7% 

4.580  x  10"25 

1.018 

S32  („,P)P32 

4.83  x  IQ'32  ±  6% 

4.785  x  10' 32 

1.009 

Ni58  (n,  p)  Co58 

1.68  x  10'32  ±  7% 

1.739  x  10“32 

0.966 

Mg24  (n,  p)  Na24 

3.25  x  10"31  ±  12% 

3.294  x  10"31 

0.987 

Al27  (n,  a)  Na114  | 

1.92  x  10'31  ±  12% 

1.849  x  10'31  j 

1.039 

Standard  Percent  Deviation  of  the  Ratios  =  2.08% 


(a)  The  resonance  (n,y)  foils  were  irradiated  inside  an  ~  40  mil  cad¬ 
mium  cover,  The  fission  foils  were  irradiated  inside  ~  40  mil 
cadmium  placed  inside  a  boron  carbide  spherical  cover  ~  1  cm 
thick.  The  other  foils  were  bare, 

(b)  All  measurements  made  at  Sandlo  Corporation  (Ref.  20).  Experi¬ 
mental  errors  listed  are  one  standard  deviation. 
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Tabli:  XXI 
SMFS  TEST  B 

SAND  II  THIRTIETH  ITERATION  RATIOS  OF  MEASURED  TO 
CALCULATED  ACTIVITIES 


Foil.  .  . 
Reaction  a 

Measured 
Activity 
(Dps/Nucleus ) 

Calculated 

Activity 

(Dps/Nucleua) 

Ratio 

Meas.  /Calc. 
Activities 

Cu63  (n,y)  Cu64 

1.63  x  10'30  ±  10% 

1.632  x  10’30 

0.999 

Au197(n,y)Au198 

2.20  x  10' 29  ±  5% 

2.204  x  10'29 

0.998 

Co59  (n,  y)  Co60 

3.41  x  10'33  ±  10% 

3.416  x  10‘33 

0.998 

Na23  (n, y)  Na24 

2.59  x  10“32  *  5% 

2.594  x  10‘32 

0.999 

In115(n,n')ln115m 

6.9?  x  10‘30  ±  10% 

7.073  x  10"30 

0.985 

Fe54  (n,p)  Mn54 

2.68  x  10‘33  ±  10% 

2.342  x  io‘33(c) 

1,  145 

Ni58  (n,p)Co58 

1.12  x  10' 32  ±  5% 

1.204  x  10“32(d) 

0.930 

S32  (n.p)P32 

3.16  x  10‘32  *  6% 

3.270  x  10"32 

0.966 

Zn64  (n,  p)  Cu64 

6.62  x  10‘31  *  10% 

6.517  x  10“31 

1.016 

Cl35  (n,a)P32 

8.95  x  10"33  ±  10% 

9.523  x  10“33 

0.940 

Al27  (n, a)  Na24 

1.48  x  10’31  ±  12% 

1.436  x  10“ 31 

1.031 

Mg24  (n>P)  Na24 

2.54  x  10'31  ±  12% 

2.517  x  10“31 

1.009 

Standard  Percent  Deviation  of  the  Ratios  =  6,09% 


(a)  The  reaonance  (n,y)  foil*  were  all  irradiated  Inside  an  20  mil  cad¬ 
mium  cover,  except  for  the  Cu  which  was  irradiated  bare.  The  other 
foils  were  bare. 

(b)  Measurements  made  at  Atomics  International  and  Sandia  Corporation, 
(Ref.  20).  Experimental  errors  listed  arc  one  standard  deviation, 
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Table  XXII 
BRGM  TEST 


SAND  II  THIRTIETH  ITERATION  RATIOS  OF  MEASURED  TO 
CALCULATED  ACTIVITIES 


Foil 

Reaction(a) 

Saturated 
Measured 
Activity 
(Dps /Nucleus) 

Saturated 

Calculated 

Activity 

(Dps/Nucleus) 

Ratio 

Meas,  / Calc. 
Activities 

1 

Au197(n.r)  Au198 

1.62  x  10"21  ±  3% 

1.621  x  10"21 

Co59  (n,  y)  Co60 

7.30  x  10“23  ±  10% 

7.304  x  10"23 

Na23  (n,  y)  Na24 

2.77  x  10"25  *  4% 

2.772  x  10'25 

0.999 

Cu63  (n,y)Cu64 

5.24  x  10‘24  ±  3% 

5.244  x  10‘24 

0.999 

In115  (n,n')In115m 

3.55  x  10"25  ±  5% 

3.494  x  10’25 

1.016 

Ni58  (n,p)Co68 

1.08  x  10-25  ±  3% 

1.170  x  10"25 

0.923 

Fa 54  (n,p)  Mn54 

9.45  x  10‘2t’  ±  5% 

9.531  x  10"26 

0.992 

S32  (n,  p)  F32 

6.22  x  10”26  *  3% 

6.291  x  10"26 

0.989 

Al27  (n,a)  Na24 

3.07  x  10’28  *  12% 

3.078  x  10"28 

0.998 

Cl35  (n,  a)  P32 

1.57  x  10'26  *  10% 

1.446  x  10'26 

1.086 

I127  <n,2n)I126 

7.62  x  10’28  ±  25% 

7.620  x  10‘28 

1.000 

Standard  Percent  Deviation  of  the  Ratios  -  3.7% 


(a)  All  foils  were  Irradiated  Inside  an  ««  20  mil  cadmium  cover. 

(b)  Measurements  made  at  IIT  Research  Institute,  (Ref.  22),  Experi¬ 
mental  errors  listed  are  quoted  limits  of  error. 
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Fe  (n,p)Mn  reaction,  whic  h  in  this  test  as  in  the  GMRW  test  (Fig.  14)  is  man¬ 
ifest  as  a  spectral  depression.  In  addition  to  errors,  low  values  of  measured 
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activity  for  the  Fe  (n,  p)Mn  reaction  can  be  attributed  to  flux  depressions  re¬ 
sulting  from  materials  (such  as  the  1/2”  of  aluminum  and  stainless  steel  for 
the  GMRW  test  and  ''*■  1/2"  Hastelloy  "C"  for  the  BRGM  test)  present  in  an  en¬ 
vironment  (see  Fig.  11).  Some  experimental  results  on  the  flux  depression  ef¬ 
fect  cf  aluminum  and  stainless  steel  are  presented  in  the  discussion  of  the  ECEL 
test  (Section  V-B-3). 

The  curves  presented  at  the  bottom  of  Fig.  16  are  three  of  those  shown  in 
Fig.  11,  and  are  included  to  help  establish  the  credibility  of  some  features  of 
the  solution  spectrum.  It  is  noted  that  the  major  depression  in  the  solution  spec¬ 
trum  in  the  vicinity  of  approximately  3  Mav  is  shifted  to  the  left  and  broadened 
compared  with  the  results  of  the  GMRW  test  (Fig.  14);  this  shift  is  just  as  might 
be  expected  to  result  from  the  presence  of  beryllium  oxide,  since  beryllium  has 
a  major  resonance  at  about  2.7  Mev. 

The  ECEL  te«t  results  will  be  discussed  in  more  detail  than  the  previous 
tests,  since  these  results  demonstrate  most  of  the  ways  in  which  the  iterative 
technique  can  be  used  to  extract  various  types  of  information  from  foil  activity 
measurements. 

The  ECEL  Core  14-13  absolute  value  of  integral  flux  as  determined  by  the 
iterative  method  was  about  a  factor  of  three  lower  than  that  estimated  on  the  basis 
of  calculations  and  a  calibrated  power  monitor;  with  lower  than  expected  foil  ac¬ 
tivities,  considerable  effort  had  to  be  expended  by  the  participating  laboratories 
to  obtain  satisfactory  activity  measurements.  Nevertheless,  the  overall  con¬ 
sistency  of  the  measurements  was  excellent.  In  general,  the  determination  of 
the  absolute  value  of  total  integral  flux  is  not  a  trivial  problem.  It  is  noted  that 
the  multiple  foil  activation  method  is  rather  unique  in  this  regard  and  can  be 
used  as  a  method  of  absolute  calibration  for  most  neutron  environments. 

For  the  WMGR.,  GMRA,  and  GMRW  tests,  there  are  definite  limits  on  the 
reliability  of  the  reactor  physics  spectral  calculations,  because  of  geometrical 
and  non-homogeneity  considerations.  Because  of  the  homogeneous  and  symmet¬ 
rical  nature  of  the  ECEL  core  region,  however,  it  would  be  expected  that  theo¬ 
retical  predictions  would  agree  better  with  measurements  for  an  experiment 
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performed  in  this  facility.  This  was  indeed  found  to  be  the  case,  for  the  relative 
integral  flux  agreement  in  the  ECEL  test  is  within  20%  over  the  entire  energy  re¬ 
gion,  compared  with  factors  of  £  or  more  in  some  of  the  previous  tests  (see  Ta¬ 
bles  XV,  XVII,  and  XIX'.  Prior  to  the  EC  EL  test,  only  in  beam  type  tests  had 
agreement  been  obtained  with  theory  to  within  20%. 

The  preliminary'  foil  activity  and  differential  flux  results  using  the  SAND  II 
code  with  12,  11,  and  9  foil  reactions  for  the  +Y  Cd-covered,  -Y  Bare,  and  +Y 
B  10-covered  locations,  respectively,  for  the  ECEL  test  are  presented  in  Ta¬ 
bles  XXIII,  XXIV,  XXV,  and  Figs.  17,  18  and  19.  Values  of  10%  and  1.96  were 
used  for  the  solution  and  foil  discard  criteria  for  all  three  SAND  II  runs.  Two- 
inch-square  voids  were  created  approximately  4  inches  from  the  core  center 
along  the  x,  y,  and  z  axes.  The  two  voids  on  the  y  axis  were  designated  as  the 
+Yand  -Y  locations.  A  description  of  the  ECEL  facility,  which  is  of  the  split- 
table  honeycomb  type,  has  been  given  by  Mountford  and  Morewitz  (Ref.  25). 

Tuttle  (Ref.  30)  (using  the  CAESAR  IV  diffusion  code  with  18-group  and  35- 
group  structures)  and  Dahl  (Ref.  32)  (using  the  S-XIII  transport  code  and  18 
groups)  made  one -dimensional  spherical  geometry  reactor  physics  calculations 
for  the  ECEL  core  14-13.  (These  calculations  did  not  distinguish  between  the 
+Y  and  -Y  locations.  )  Dahl  obtained  two  results  using  different  sets  of  cross 
section  evaluations.  Tuttle's  and  Dahl's  results  are  given  in  Table  XXVI,  Tuttle's 
18-group  results  were  used  to  generate  an  input  spectral  approximation  for  the 
SAND  II  runs  shown  by  the  dashed  curves  in  Figs.  17  and  19. 

Tuttle  (Ref.  30)  also  used  the  AILMOE  code  to  calculate  a  spectrum  with  de¬ 
tailed  structure  for  core  14-13,  It  was  suggested  that  the  results  of  these  calcu¬ 
lations  be  used  only  to  indicate  where  some  detailed  structure  might  appear  in 

-3  t 

the  spectrum  between  about  10  and  3.6  Mev,  and  that  the  previous  CAESAR  IV 
(either  the  18-group  or  subsequent  35-group)^  results  be  relied  upon  for  the  cor¬ 
rect  overall  distribution  of  the  neutron  flux.  The  AILMOE  code  results  are 


>!<At  the  time  that  this  report  was  prepared,  not  all  of  the  experimental  foil  acti¬ 
vation  data  had  been  completely  analyzed. 

■fBelow  10*3  Mev  and  above  3.6  Mev,  the  AILMOE  code  assumes  an  E  form  and 
a  fission  form,  respectively. 

SDifferent  sets  of  cross  sections  were  used  for  the  two  calculations. 


-78- 


Table  XXIII 

SAND  II  SOLUTION  RESULTS  OBTAINED  AFTER  20  ITERATIONS  FOR  THE  ECEL  TEST  +Y 

Cd-COVERED  LOCATION  ’ 


d X 

.2  u 

3  3 

S  rt 
4)  «> 

QS 


xl 

0  d  > 

u  (j  U 

H  U 

rt  < 

U 


O  ro 

o  so 


CO 

in 


o 

I 


0 

i-H 

r- 

0 

r- 

rO 

in 

tfi 

CO 

"H 

fM 

NO 

in 

M 

pH 

cn 

d 

1 

d 

1 

d 

d 

d 

<N 

n J 

1 

d 

1 

in 

m 

ph 

0 

*H 

B 

o 

2 


>- 

■a  « 
y  ■ “ 
«  E 


v**i 

a  is  i 

Jp4 


0) 

a 

X) 


d 

o 

3  3 
o  u 
(Z  m 

^  jj 

ftj 


x 

vO  *o 

CO  vO  r-H  00 

^  O  H  H 

odd  d  *-i 


^  in  Tf  in  in  h  o 

*••••* 
in  o  r^.  r-  Tf 


in 

*© 


sO  sO 

■o  'o 


"t  Tf 

'o  'o 


IT) 

‘o 


•51 3  f! 

pi! 

rt  rt  ^  ia 
MU  ft 

X 


X 

CM 

CM*  Tfi 


X  X 

oo  oo 
■ 


X 

N 


X  X 

vO  o 
*  • 

in  in 


rt<  r«. 


co 

oi  o 


-H  P"*  P"!  pH  pH^CphHhHpHpHpHfH 


o  'o  'o  'o  'o  'o  'o  '©  'o  '©  1 


o  © 


X  X 


sO 

r~ 


X 

it 


CO 
vd 
O'  CO 


xxxxxxxxx 

minN'Ot'-'-HrOcOi-n 


CM 

"M  m 
in  oo 


•o 

CO 


t-p 

IM 

vO 


N  fO  M 


in  m 


in 

■ 

r- 


6^ 


fiPo 

[Mphph[pcoco(ooocoip..-h,-< 


«  *H  -H  -H  -H  -H 
so  in  m  in  vO  in 


■H  -H  ■«  +)  ■«  -H 

-  n-  n  oo  o 


o  'o  'o  'o  'o  'o  'o  '©  'o  'o  'o  'c 


XXX 


cm  in 
cn 

0-  -H 


O' 

sC 

© 


X 


XXXXXXXX 


nj  O 

CM  CO 

in  o 


o 

N 

•O 


© 

00 

N 


rt<  o 
f».  so 
in  cm 


00 

pjl 

CO 


n  n  n 


S© 


^  & 
pi  x  ""fl 
X  U  <3  x  x  X 

1}  oo  oo  O  0  U 

£  £  Z  °  * 

s  ^  <  U  U 

N  c? 


x  x  x  x  x  x 

o  o  u  u  u  u 


.  oo  rt*  © 
d  m  in  .  m  cm 

&i  u  : 


<M  nj 

*o<  “h 


g  £  3,  a  g 

n-  fp  w  w 

oo  o  on  ^  rn  m  oo 

rO  '“H  -  - 

<m  d 

D  <  <J  0  U 


a  a  3 

*•»  *w  _ 


ft  a  ©  cm 


CM 

CO 


cm  ®  in  p,  n- 

B  0  3  cm  O)  -B  lH  ,k  Io  23 


3  D  H  2  h 


M 


(A 

■d 

■? 

t 

X 

0) 

M 

3 

U) 

rt 

IU 

2 

PH 

O 

fl 

0 

td 

#5! 

<s 

X 

h 

rt 

x 

£ 

U3 


4) 

pH 

»“4 

HH  4) 

W  a. 

iH 

^  rt 
oo 

-H  -0 

(V 

•  4-J 
«  2 
ej  ^ 

—  U) 
*H 

0  ^ 
*j  Ih 
rt  V 
Ih  ^ 
0  rt 

hQ  pH 

•  rt  d 

f?  J 
2h  E 
u  rt  'r! 
#  fl  J 

4J  O  D 

4)  3  a 
X  «  X 

U 

i  b- 

X  S© 
C  O 
rt  W)  , 
w  u  m* 
<q  « 

4>  nT 

»  »H 

rt  rt  rt  r 

u  y  G  o 

0  'H 
»  «,h" 
>H  H  H  (w 

Jq  Jc  rt  x 
«  **  rt  o 

u  a. 
<u  u 

«  “3o 


fl  c 


1  a 

o  m 

*  d 


«  rt 

|! 

q  to 


is  ,S  1 <  ’2 

>  ?  +j  S 

8  «  rt 

a  0  X  _  , 

a  V  S 

■g  x  S  , .« 

6  ’o  «>  'S  rt 

«  m  a  ctJ  S 

5  rt 

_  C3  w  ^3 

*j  *j  a)  ^  x 

W  W  H  St 

Ih  y  3  Jp  J 

4>  4)  S2  *3  Y 

Ih  i*h  ™  li  ^ 

S 


rt  -o  u 


-79- 


TJ 

J8  >> 

6  3's 

O  3  ,> 


m  v  'f  eo 
M  HI  h  M 


vO 

o> 

<M 

on 

CO 

6? 

vo 

Tf 

r- 

00 

CO 

D 

o 

00 

d 

rsJ 

o} 

vO 

d 

rO 

X  y 

^  Tp  ”  sO 

(SJO00~<rtOT}<T}'if|^00' 

oo<— |O'-Hir)s0r^*r^r^-T^ 


Lfi  sO  ^  iJH  i/) 

*o  *©  'o  'o  '© 

"H  —I  — I  — < 

X  X  X  X  X 

*m  ^  «  in  in  i-J  i-4  —<  — I  -o" 


’o 

•o 

d 

•o 

•o 

‘o 

d 

•o 

d 

'© 

•o 

pH 

p-H 

-H 

-h 

H 

<“H 

pH 

pH 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

O 

O 

CO 

<D 

pH 

no 

ON 

r^ 

OJ 

*NJ 

N 

O 

pH 

h- 

f\) 

CO 

ON 

o 

S o 

N 

N 

LO 

<M 

ID 

D 

pH 

Tt* 

CO 

d 

d 

pH 

— < 

d 

pH 

-H 

pH 

• 

h- 

-I 

-h-h-h-h-h-h-h-h-h** 


'o 

o 

d 

d 

'o 

d 

'o 

'© 

d 

'o 

d 

pH 

*-H 

pH 

►H 

<-H 

pH 

pH 

pH 

pH 

pH 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

vO 

CO 

OJ 

ON 

o 

O 

o 

sO 

CO 

On 

CO 

fN- 

OJ 

•o 

OJ 

s0 

H 

h- 

CO 

pH 

0s 

CO 

pH 

if) 

PH 

sO 

o 

D 

r0 

ON 

co 

CO 

CO 

pH 

d 

d 

— ■ 

pH 

d 

H 

5  tf  o:  s 

6  ^  ^  ^  h  h  hi  U  |  ^  « 

*  ^  5  *  ^  h  \  *  *?  a  ^ 

rT  C3  C  C  g'fl' 


^  rvj 


[N* 

00  ON 

ON 

CO 

D 

00 

<n 

CO 

CO  pH 

D 

sO 

cO 

rO 

OJ 

n  d 

0 

d 

CvJ 

cvj 

pd 

< 

U 

U 

D 

Id 

fH 

’  !!j  *-■  t 

G  ~  s 


*i  v  i> 
*  a!  « 

a;  «—  xl 

x  w  Tj 

m  J: 

W  C 
£  £  rt 

C  In  4J 

(DOM 

£2  D 

2  «  S 


C  ^  4i 
D  ,,  C 
M  ®  O 

3  =5  _. 

W  D  X 
(Tj  -W  <U 
J,,  *J  H 

i}  rt  'A 

S  CQ  a 


SAND  II  SOLUTION  RESULTS  OBTAINED  AFTER  18  ITERATIONS  FOR  THE  ECEL  TEST 

+Y  BORON  10 -COVERED  LOCATION 
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NEUTRON  DIFFERENTIAL  FLUX  I  W-sec-Mevi 
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-  -  -  -  18-GROUP  CAESAR  IV  CALCULATION. 

-  SAND  II  RESULT, 

.  L  ENERGY  LIMITS  WITHIN  WHICH  90%  OF  THE 
ACTIVITY  IS  PRODUCED  FOR  THE  INDICATED 
FOIL. 

-THIS  Au  FOIL  HAD  AN  0.5  mil  Au  COVER. 


M 


NEUTRON  ENERGY (Mev) 

Figure  17.  SAND  II  Solution  Spectrum  After  20  Iterations  for  EC  EL 
+  Y  Cd-Covered  Location 
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-  -  -  16-GROUP  CAESAR  IV  CALCULATION. 

—  SAND  II  RESULT. 

-I  I*-  ENERGY  LIMITS  WITHIN  WHICH  90%  OF  THE 
ACTIVITY  IS  PRODUCED  FOR  THE  INDICATED 
FOIL, 
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Figure  19.  SAND  II  Solution  Spectrum  After  18  Iterations  for  ECEL 
Test,  +Y  Boron  10-Covered  Location 


Table  XXVI 

ECEL  CORE  14-13,18-GROUP  CALCULATIONS 
(NORMALIZED  INTEGRAL  FLUX  ABOVE  THE  SPECIFIED  ENERGY) 


Energy  (Mev) 

CAESAR  IV<a> 

S-XIII<b) 

S-XHI(c> 

Ratio(a) (b) (c) (d) 

IV/ XIII 

3.0 

0.02534 

0.02260 

0,02552 

1. 121 

1.4 

0. 1130 

0.09665 

0.09289 

1.  169 

0.9 

0.  1688 

0.1448 

0.1444 

1. 166 

0.4 

0.3029 

0.2607 

0.2  34  3 

1.  162 

0.1 

0.5282 

0.4628 

0.4101 

1.  141 

1.70xl0"2 

0.7344 

0.6715 

0.6926 

1.094 

3. 35x10"  3 

0.8779 

0.8379 

0.8414 

1.048 

4.54xl0'4 

0.9657 

0.9548 

0.9571 

1.011 

6. 14x10*  5 

0.9955 

0.9956 

0.9934 

1.000 

2. 26x10"  5 

0.9989 

0.99935 

0.99933 

0.999 

8.315xl0“6 

0.99960 

0.999918 

0.999895 

1.000 

3.059xl0“6 

0.999764 

0.999983 

0.999930 

1.000 

1.125xl0-6 

0,999893 

0.999989 

0.999959 

1.000 

4. 14x10*  7 

0.9999728 

0.9999956 

0.999989 

1.000 

1. 52xl0"7 

0.9999952 

0.9999979 

0.9999988 

1.000 

5.60x10'® 

0.99999955 

0.9999988 

0.9999989 

1.000 

2. 52x1 0*8 

0.99999996 

0.9999994 

0.9999995 

1.000 

9.26xl0*9 

1.00000000 

1.0000000 

1.0000000 

1.000 

(a)  Calculated  at  Atomics  International  (Ref.  30),  Jan. ,  1967 
(Diffusion  code  results). 

(b)  Calculated  at  Battelle  Northwest  (Ref.  32),  April,  1967 
(Transport  code  results,  Program  S-X1II  using  cross  sections 
from  reference  34.  ) 

(c)  Same  as  (b),  but  using  the  cross  sections  for  Hanford  revised 
Gam  (Ref.  32). 

(d)  Ratio  of  the  2nd  and  3rd  columns. 


given  in  Fig.  20.  The  structure  shown  in  Fig.  20  is  not  evident  in  the  SAND  II 
solutions,  Figs,  17,  18  or  19. 

The  SAND  II  solution  values  of  absolute  integral  flux  (4*  >  10"^  Mev)  for  the 

q  q 

+  Y  Cd-cuvered,  -  Y  Bare,  and  +  Y  B  10 -cove  red  locations  are  1.17x10,  1,12x10, 

9  2 

and  1.18  x  10  n/cm  -Sec  respectively,  at  a  power  of  150  watts.  (As  previously 
indicated,  these  results  are  based  on  sets  of  12,  11,  and  9  foil  reactions  re- 

9 

spectively.  )  An  initial  estimate  based  on  other  data  yielded  a  value  of  ■—  3  x  10 
n/cm^-sec.  The  theoretical  and  measured  normalized  integral  flux  results  at 
energy  points  used  for  the  35-group  calculation  are  presented  in  Table  XXVII, 

As  previously  stated,  the  agreement  between  measurement  and  theory  is  within 
~  ±  20%  at  all  energy  points. 

For  the  +Y  B  10-covered  test,  the  boron  10  carbide  sphere  would  be  ex¬ 
pected  to  cause  a  flux  depression  that  would  reduce  the  ECEL  core  fission  rate 
in  its  immediate  vicinity.  The  integral  flux  results  presented  in  Table  XXVII  do 
show  such  an  effect.  The  SAND  II  code  calculates  the  attenuation  effect  of  the 
boron  10  carbide  sphere  only  approximately  (and  results  would  therefore  not  be 
expected  to  be  extremely  accurate),  because  of  approximations  required  re¬ 
garding  geometry,  scattering  equilibrium,  etc.  Nevertheless,  the  -Y  B  10- 
covered  results  appear  to  be  very  consistent  with  the  +Y  Cd-covered  and  -Y 
bare  results. 

In  Table  XXVII,  Tuttle's  35-group  CAESAR  IV  results  are  used  as  the  theo¬ 
retical  spectrum.  In  addition,  the  SAND  II  results  are  given  for  seven  energy 
points  that  are  not  included  in  the  35-group  results,  but  were  used  for  the  18- 
group  calculations  of  Tuttle  and  Dahl.  The  SAND  II  results  at  those  seven  points 
may  therefore  be  compared  with  the  corresponding  values  from  the  18 -group 
calculations,  given  in  Table  XXVI.  It  is  noted  that  these  latter  calculations  are 
for  a  radius  of  0  cm  (i.  e.  ,  the  core  center),  whereas  the  35-group  calculation 
is  for  a  radius  of  10  cm,  which  corresponds  to  the  actual  location  of  the  foil 
sets;  however,  because  of  the  homogeneous  nature  of  the  core  region,  the  nor¬ 
malized  results  should  not  be  significantly  different. 

In  selection  of  individual  foils  for  flux  measurements  in  a  given  environment, 

the  energy  sensitivity  limits  in  the  anticipated  spectrum  are  very  important. 

The  SAND  II  results,  Tab.es  XXIII,  XXIV,  and  XXV,  and  Figs.  17  and  19  illus  - 

59 

trate  this  fact  well;  for  the  ECEL  solution  spectrum  in  Cd  (or  bare),  the  Co  (n,y) 
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Table  XXVII 


ECEL  CORE  14-13  INTERLABORATORY  REACTOR  TEST  INTEGRAL 

FLUX  COMPARISON 

(SAND  II  Zcf'1,  7t*\  and  18^'  Iteration  Normalized  Integral  Flux 
Comparison  For  the  +Y  Cadmium-Covered,  -Y  Bare,  and 
+  Y  Boron  10-Covered  Locations  Respectively) 

(Page  1  of  2) 


Energy  (Mev) 

Theoretical  (a) 

Experimental 

Ratio  (b) 
Exp.  /Theo. 

+  Y  Cd 

-Y  Bare 

+  Y  B  10 

9.26xl0'9 

1.0000000 

1.0000 

1.0000 

1,0000 

1.000 

4. 14x10"  7 

0.9999828 

1.0000 

1.0000 

1.0000 

1.000 

6.83xl0"7 

0.999876 

1.0000 

1.0000 

1.0000 

1.000 

1.13xl0‘6 

0.999724 

0.9999 

0.9999 

0.9997 

1.000 

1.86 

0.999527 

0.9998 

0.9998 

0.9996 

1.000 

3.06 

0.99935 

0.9997 

0.9996 

0.9995 

1.000 

5.04 

0.99924 

0.9996 

0.9994 

0.9995 

1.000 

8.32 

0.99915 

0.9995 

0.9992 

0.9994 

1.000 

1, 37xl0'5 

0.9989 

0.9992 

0.9990 

0.9985 

1.000 

2.26 

0.9986 

0.9986 

0.9985 

0.9965 

1.000 

3.73 

0.9977 

0.9978 

0.9973 

0.9930 

1.000 

6.14 

- 

0.9944 

0.9939 

0.9895 

- 

8. 14 

0.9956 

0.9928 

0.9914 

0.9874 

0.997 

l.OlxlO'4 

0.9915 

0.990 

0.989 

0.983 

0.998 

1.67 

0.9849 

0.982  j 

0.980 

0.957 

0.997 

2.75 

0,9745 

0.974 

0.970 

0.947 

0.999 

4.54 

0.9593 

0.959 

0.950 

0.941 

0.9999 

7.49 

0.9396 

0.943 

0.928 

0.934 

1.004 

1.23xl0'3 

0.9148 

0.922 

0.900 

0.926 

1.008 

2.04 

0.8838 

0.897 

0.867 

0.918 

1.015 

3.35 

0.8468 

0.869 

0.831 

0.894 

1.026 

5.53 

0.8050 

0.833 

0.794 

0.828 

1.035 

9.12 

0.7662 

0.801 

0.757 

0.794 

1.045 

1.50xl0"2 

0.7152 

0.761 

0.717 

0.768 

1.064 

1.70 

- 

0.751 

0,706 

0.760 

- 
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Table  XXVII 


ECEL  CORE  14-13  INTERLABORATORY  REACTOR  TEST  INTEGRAL 

FLUX  COMPARISON 

(SAND  II  20^,  7^\  and  18^  Iteration  Normalized  Integral  Flux 
Comparison  For  the  +Y  Cadmium-Covered,  -Y  Bare,  and 
+Y  Boron  10-Covered  Locations  Respectively) 

(Page  2  of  2) 


Energy  (Mev) 

Theoretical  (a) 

Experimental 

Ratio  (b) 
Exp.  /Theo. 

+Y  Cd 

2.48 

0.6589 

0.717 

m 

4.09 

0.6049 

0.669 

0.620 

mm 

6.74 

0.5399 

0.612 

0.564 

□H 

1.00x10" 1 

- 

0.5574 

0.512 

0.587 

- 

l.llxlO"1 

0.4755 

0.543 

0.500 

0.570 

1.142 

1.83 

0.4073 

0.459 

0.422 

0.467 

1.127 

3.02 

0.3311 

0. 364 

0.337 

0.338 

1.099 

4.00 

- 

0.309 

0.292 

0.263 

- 

4.98 

0.2514 

0.267 

0.256 

0.213 

1.062 

8.21 

0.1724 

0.187 

0.186 

0.138 

1.085 

9.00 

- 

0.176 

0.175 

0.128 

- 

1.35x10° 

0.1031 

0.130 

0.125 

0.0900 

1.261 

1,40 

- 

0.123 

0, 120 

0.0867 

- 

2.23 

0.04545 

0.0393 

0.0422 

0.0314 

0.865 

3.00 

- 

0.0197 

0,0213 

0.0166 

- 

3.78 

0.01504 

0,0160 

0.0126 

0.0102 

1.064 

6.07 

0.002737 

0.00244 

0.00273 

0.00219 

0.891 

l.OOxlO1 

0 

0.000206 

0.000209 

0.000108 

- 

a)  35-Group  CAESAR  IV  Calculation  (Ref.  30). 

These  results  are  for  a  radial  distance  of  10  cm,  which  corresponds  to 
the  actual  location  of  the  foil  sets  for  the  X,  Y,  and  Z  locations. 

b)  The  ratios  given  are  for  the  +Y  Cd  location.  The  ratios  for  the  -Y  bare 
location  are  closer  to  unity,  in  general,  than  for  the  +Y  Cd  location. 
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197 

reaction  behaves  in  the  same  manner  as  does  Au  (n,y)  in  a  1/E  environment; 
i.  e.  ,  90%  of  the  activity  is  produced  by  a  very  narrow  energy  region  around  the 
major  resonance.  This  behavior  makes  it  clear  that  resonance  detectors  can¬ 
not  be  used  indiscriminately  as  flux  monitors,  without  knowledge  of  the  energy 
sensitivity  range  of  the  detector  in  question,  in  the  particular  environment  in 
which  it  is  to  be  used. 

It  is  instructive  to  examine  some  additional  observations  regarding  the 
dips  and  peaks  in  the  three  ECEL  solution  spectra: 

a)  For  the  SAND  II  data  reduction,  the  fission  foil  results  of  Armani 
(Ref.  18)  for  the  +Y  Cd-covered  and  -Y  bare  locations,  and  those  of  Tochilin 
(Ref.  19)  for  the  +Y  B  10-covered  locations,  were  included  as  input  data. 
Armani's  results  are  based  on  counting  of  reaction  products,  while  Tochilin's 
are  based  on  counting  of  fission  tracks.  The  high  value  of  the  integral  spectrum 
ratio  (Table  XXVII)  at  1.35  Mev  and  the  spectral  peaks  indicated  by  the  arrows 
in  Figs.  17  and  18  at  that  energy,  appear  to  be  due  to  the  Th  (n,f)  reaction, 
which  has  high  ratios  of  measured  to  calculated  activity  for  both  the  +Y  Cd- 
covered  and  -Y  bare  locations  (see  Tables  XXIII  and  XXIV).  , 

This  same  type  of  behavior  was  observed  in  the  previous  GMRW  test, 
232 

as  seen  in  Fig.  1.  The  Th  (n, f)  results  for  this  test  also  were  provided  by 

232 

Armani  by  radioactivity  counting.  When  the  Th  (n,  f)  reaction  was  omitted 

from  the  GMRW  test  results,  the  SAND  II  solution  flux  spectrum  changed,  the 

peak  was  greatly  reduced,  and  the  spectral  dip  on  the  low  energy  side  of  the 

original  peak  shifted  up  in  energy,  as  shown  in  Fig.  21.  The  peak  observed  in 

each  of  the  spectra  of  Figs.  17  and  18  (shown  by  the  arrows)  does  not  appear  in 

Fig.  19  at  the  same  energy,  but  is  shifted  up  in  energy  as  in  Fig.  21,  (It  is 

232 

again  noted  that  the  Th  (n,  f )  results  for  the  +Y  B  10  test  (Fig.  19)  were  pro¬ 
vided  by  Tochilin  by  counting  fission  tracks.  ) 


❖Regarding  the  behavior  of  the  gold  detector  in  a  spectrum  such  as  the  ECEL,  it 
might  be  noted  that  in  the  +Y  Cd-covered  location,  a  gold  covered  gold  (sand¬ 
wich)  detector  was  used  (Table  XXIII),  but  the  relative  hardness  of  this  envi¬ 
ronment  caused  the  gold-covered,  cadmium-covered,  and  bare  gold  detectors 
to  have  essentially  the  same  energy  response  range.  Since  the  flux  at  and  be¬ 
low  the  major  gold  resonance  (~  5x  10"°  Mev)  is  relatively  small  compared  to 
a  1/E  spectrum,  the  cadmium  cover  (which  attenuates  the  thermal  neutrons 
preferentially)  and  the  gold  cover  (which  attenuates  preferentially  in  the  region 
of  the  gold  resonance)  had  relatively  little  effect. 
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These  observations,  as  well  as  subsequent  discussion  with  the  partic- 

232 

ipants,  suggest  chat  the  i'h  '(n,  f)  fission  rates  as  determined  by  the  counting 

of  reaction  products  may  be  excessively  high  compared  with  those  determined 

by  fission  track  counting,  probably  due  to  the  uncertainty  in  the  fission  yield  of 

99 

the  Mo  reaction  product  that  is  counted.  The  possibility  of  cross  section  er¬ 
ror,  however,  cannot  entirely  be  ruled  out. 

b)  Examination  of  the  input  spectrum  used,  shown  in  Figs.  17  and  19  as 
the  dashed  curve,  may  raise  conjecture  that  the  peak  discussed  above  is  pos¬ 
sibly  caused  by  the  fact  that  the  input  spectrum  includes  a  small  peak  in  that 
energy  region.  This  is  most  definitely  not  the  case,  however.  Arepeat  SAND  II 
run  was  made  for  the  +Y  Cd-covered  location,  using  a  smooth  input  form  with 
no  structure  in  this  region;  the  peak  produced  in  the  solution  for  this  run  was 
essentially  the  same  as  that  shown  in  Fig.  17  in  magnitude  and  only  very  slightly 
different  in  shape.  The  results  of  this  comparison  are  shown  on  an  expanded 
scale  in  Fig.  22.  It  can  be  seen  that  the  difference  in  shape  consists  of  the  pres¬ 
ence  or  absence  of  the  remnant  of  the  small  input  peak  at  about  1.7  Mev,  but  that 
the  generation  of  the  solution  peak  {at  1.35  Mev)  is  independent  of  the  input  peak. 

Since  the  activity  results  for  the  Na  °(n,y),  In  (n,n')>  In  (n,y), 

56  27  24 

Fe  (n,  p),  A1  (n,a)  and  Mg  (n,  p)  were  not  available  for  the  current  computer 
runs,  the  energy  region  from  I  to  8  Mev  did  not  have  proper  foil  coverage. 
Additional  structure  in  the  solution  differential  spectra  in  this  energy  region 
may  be  detected  when  these  reactions  are  used.  There  is  no  doubt  as  to  whether 
such  structure  could  actually  exist.  The  plausibility  of  the  SAND  II  solution  for 
the  GMRW  test,  shown  in  Fig.  21,  is  strongly  supported  by  the  representative 
spectra  shown  in  Fig.  11  (which  show  the  effects  different  materials  actually 
do  have  in  causing  spectral  depressions).  Further,  there  is  no  question  about 
the  ability  of  the  SAND  II  code  to  detect  such  structure  if  accurate  c  ros  s  sections, 
accurate  activities,  and  sufficient  foil  coverage  are  obtained.  The  current  prob¬ 
lem  is  rather,  with  existing- uncertainties,  to  properly  analyze  SAND  II  solution 
spectra  to  determine  whether  the  structure  produced  is  real  or  is  the  result  of 
errors  in  cross  sections,  measured  activities,  foil  self-shielding  effects,  or 
inadequate  foil  coverage  of  the  entire  energy  range. 
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BASED  ON  12  FOIL  REACTIONS  AS  INPUT  -  ADDITIONAL 
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Al27  (n,  p),  AND  Mg24(n,  p)  RESULTS  ARE  USED. 
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c)  In  Fig.  17,  two  small  dips  in  the  10  Mev  decade  (one  between  4  and 

5  9 

5  and  one  between  7  and  8)  very  closely  coincide  with  two  Co  (n,  y)  resonances. 

(Figure  19  shows  the  artificial  peaks  produced  by  these  resonances  in  the  +Y  B 

10-covered  test  very  clearly.  )  It  might  appear  on  first  examination,  therefore, 

59 

that  these  two  dips  are  caused  by  the  Co  (n,y)  reaction  --  i.  e.  ,  by  a  negative 
error  in  the  measured  activity  for  that  foil.  This  is  implausible,  however,  for 

two  reasons.  Firstly,  the  flux  at  the  major  Co  (n,  y)  resonance  (between  1  and 

-4 

2x10  Mev)  shows  no  such  behavior,  whereas  any  error  in  the  measurement 

would  be  propagated  most  strongly  at  the  major  resonance.  Secondly,  the  90% 

59 

activity  limits  plotted  in  Fig.  17  show  that  the  Co  (n,y)  reaction  has  negligible 

-  3 

sensitivity  in  the  10  Mev  decade,  and  could  therefore  not  be  expected  to  have 
any  effect  on  the  solution  spectrum  in  that  region. 

Further  examination  shows  that  the  dips  in  question  are  only  two  of 
63 

five  that  coincide  with  Cu  (n,y)  resonances.  The  fact  that  all  five  dips  are 

quantitatively  similar  strongly  suggests  that  these  dips  are  indeed  caused  by  a 

6  3 

somewhat  low  activity  for  the  Cu  tn,y)  foil  and  are  not  real.  This  is  another 
example  of  how  the  generation  of  structure  in  the  solution  differential  spectrum 
can  be  used  to  advantage  to  evaluate  the  reliability  of  foil  activity  measurements. 

238 

d)  In  discussions  with  Armani  (Ref.  18),  it  was  learned  that  the  U 

2 

foil  used  was  ~  240  mg/cm  thick.  In  the  GMRW  test,  Armani  had  experimen¬ 
tally  determined  that,  in  a  1/E  environment,  the  self-absorption  correction 

factor  for  such  a  foil  is  3.3.  Therefore,  the  measured  activity  for  a  240  mg/cm2 
238 

U  (n,y  )  foil,  without  self-absorption  correction,  will  be  substantially  too  low 

to  be  consistent  with  the  other  foils  used  (although  not  necessarily  by  a  factor  of 

3.3  in  the  ECEL  test,  since  it  is  not  a  1/E  environment).  In  striking  agreement 

with  this  prediction  is  the  observation  that  nearly  all  the  dips  in  the  solution 

spectrum  between  6  x  10  ^  Mev  and  3  x  10  ^  Mev  in  both  Fig.  17  and  Fig.  18 

238 

correspond  to  resonances  for  the  U  (n ,y  )  reaction.  (This  reaction  was  not 

included  in  the  +Y  B  10-covered  test  of  Fig.  19.)  In  an  environment  such  as  the 
238 

ECEL  (in  which  U  is  uniformly  distributed),  it  is  to  be  expected  that  some 

2  38 

portion  of  all  these  dips  is  real,  caused  by  the  resonances  in  the  U  present 
in  the  environment. 

A  very  interesting  hypothesis  emerges  from  this  observation:  if  the 
2  38 

measured  activity  for  the  U  (n  ,Y )  reaction  is  carefully  determined  for  a  given 
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Z  38 

thickness  of  foil  in  any  environment  where  U  is  not  in  the  immediate  vicinity 

of  the  foils,  it  is  possible  (starting  with  the  measured  activity)  to  simply  increase 

the  SAND  II  input  activity  for  that  reaction  until  all  the  dips  due  to  its  resonance 

disappear  in  the  solution  spectrum;  the  ratio  of  input  activity  to  measured  activity 

needed  to  accomplish  this  would  then  be  the  self-shielding  correction  factor  for 
238 

U  (n  ,y)  for  the  environment  under  consideration.  Obviously,  for  an  environ  - 

238 

ment  where  U  is  present,  only  an  upper  limit  for  the  self-shielding  correction 

factor  can  be  obtained  in  this  manner.  (This  technique  was  used  to  determine  a 

238 

self-shielding  correction  factor  upper  limit  of  ~  1,45  for  the  U  (n ,y)  foil  in  the 
+Y  cadmium-covered  location.  ) 

This  method  of  studying  and/or  determining  eelf-shielding  factors 
for  a  given  environment  can,  at  least  in  principle,  be  extended  to  any  resonance 
detector  reactions,  thereby  greatly  reducing  a  very  important  source  of  uncer¬ 
tainty  in  thick  resonance  foil  activiation  techniques.  One  obvious  requirement 
for  this  method  is  that  enough  other  foils  be  used  which  are  largely  redundant 
in  energy  range  of  sensitivity  with  the  foil  under  investigation  so  that,  determina¬ 
tion  of  the  true  spectral  shape  is  not  strongly  dependent  on  that  foil.  (Whether 
or  not  this  latter  condition  is  satisfied  can  be  tested  by  making  a  SAND  II  run 
with  the  foil  in  question  removed  --  except  for  the  dips  corresponding  to  the  res¬ 
onance  structure  of  the  foil  being  studied,  the  solution  spectrum  should  be  the 
same  with  and  without  the  foil.  ) 

e)  There  is  a  large  peak  in  the  Fig.  17  solution  spectrum  after  20  itera¬ 
tions  (and  only  an  incipient  one  after  7  iterations  in  Fig.  18)  that  coincides  in 

197 

energy  with  the  main  Au  (n, Y)  resonance,  In  this  case,  the  appearance  of  a 
peak  would  tend  to  suggest  a  positive  error  in  the  measured  activity,  since  the 
energy  sensitivity  region  for  this  reaction  does  extend  over  the  peak  for  this  so¬ 
lution  spectrum.'*  On  the  other  hand,  in  Fig.  19,  the  large  dip  in  the  solution 

197 

spectrum  at  this  same  Au  (n,7)  resonance  energy  is  of  negligible  significance 
with  respect  to  the  measured  activity,  since  the  energy  sensitivity  region  for 

Au  (n,y)  in  that  spectrum  does  not  extend  below  ~  2  x  10  Mev.  In  fact,  the 

-4 

entire  solution  spectrum  below  ~  1  x  10  Mev  in  Fig.  19  should  be  considered 

*The  peak  may,  however,  be  due  to  a  problem  in  the  SAND  II  code's  approxi¬ 
mate  calculation  of  the  attenuation  for  the  gold-covered  gold  sandwich  detector, 
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to  be  largely  meaningless  since,  for  this  particular  spectrum  in  the  region  of 
the  boron  10  carbide  sphere,  none  of  the  reactions  used  were  significantly  sen¬ 
sitive  below  that  energy.  It  is  thus  seen  to  be  very  important  to  consider  the  re¬ 
action  energy  sensitivity  limits  (for  the  solution  spectrum)  tabulated  by  the 
SAND  II  code  in  evaluating  solution  spectral  structure. 

f)  Additional  structure  was  introduced  into  the  ECEL  solution  flux  spec- 

35 

tra  for  the  +Y  Cd-covered  and  the  -Y  bare  locations  by  the  Cl  (n,a)  and 

Zn84(n,  p)  reactions;  however,  the  SAND  II  code  discarded  both  reactions  for 

both  locations  because  of  excessive  deviation  (>1.96  standard  deviations  high) 

from  consistency  with  the  solution  spectrum  (results  for  these  reactions  were 

not  yet  available  for  these  runs  for  the  +Y  B  10-covered  test).  Subsequent  checks 
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showed  the  presence  of  an  impurity  in  the  Cl  samples,  which  justified  the  dis¬ 
carding  of  this  reaction.  Because  of  the  very  low  activity  of  the  Zn  samples  and 
problems  associated  with  measuring  the  Cu^4  activity  (see  Section  V-A-2),  the 
Zn84(n,  p)Cu^*4  results  were  not  considered  to  be  very  accurate,  which  justified 
the  discarding  of  this  reaction  also. 

3.  Flux  Depression  Effects 

Because  of  a  suspected  fluxdepression  in  th£  GMRW  test  (due  to  the  pres¬ 
ence  of  l/4"  thicknesses  of  aluminum  and  stainless  steel,  as  well  as  water,  around 
the  foil  set  used),  a  test  was  performed  in  the  ECEL  minus  Z  location  with 
Fe54(n,p)Mn54,  Ni58(n,  p)Co58,  and  Fe56(n,  p)Mn56  foil  detectors  sandwiched  be¬ 
tween  two  1/2"  total  thicknesses  of  aluminum  and  stainless  steel  (1/4"  x  2  1  x  2 
pieces  each  of  A1  and  SS)  to  ascertain  the  degree  of  flux  depression.  Measured 
activities  for  these  sandwiched  foils  relative  to  foils  located  at  the  center  of  the 
outside  surfaces  of  the  1/2"  thick  aluminum  and  stainless  steel  sandwich,  were 
reduced  by  8.4  (±6)%  for  Fe^(n,  p),  7.8  (±3)%  for  Fe^4(n,  p),  and  7.1  (±3)%  for 

CO 

Ni  (n,  p).  Relative  to  the  average  of  activities  for  these  reactions  previously 
measured  in  the  center  of  2"  x  2"  x  2"  voids  (with  no  aluminum  and  stainless 
steel  sandwich),  the  measured  reductions  were  25.6  (±15)%,  10.5  (±5)%,  and 

12.6  (±3)%  respectively.  The  more  pronounced  effect  for  the  Fe^(n,  p)'  may  be 

________  5  A 

>!<Part  of  this  larger  difference  for  Fe  (n,  p)  can  be  attributed  to  the  presence 
of  manganese  impurity  in  the  iron  foil  (see  Section  V-A-2). 
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consistent  with  the  fact  that  it  is  the  most  abundant  of  these  three  isotopes  in 
stainless  steel.  That  is,  the  presence  oi  a  flux  depressant  material  might  be 
expected  to  have  the  strongest  effect  on  a  detector  of  the  same  material. 

Davey  and  Amundson  have  made  somewhat  analogous  measurements  for 

~  1/4"  thick  spherical  shells  of  graphite,  sodium,  aluminum,  iron,  stainless 

steel,  lead,  and  depleted  Uranium  aroundfission  chambers  {Ref.  35).  They  report 

changes  in  fission  count  rate  ratios,  on  adding  the  different  shells,  of  ~  1.0  to 

6.6%;  the  low  and  high  values  correspond  to  sodium  and  depleted  uranium  shells 

respectively.  The  high  value  for  the  uranium  is  again  consistent  with  the  ex- 
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pectation  that  the  fission  chamber  (back-to-back  halves  coated  with  U  and 
235 

U  ,  respectively)  will  be  most  strongly  affected  by  the  presence  of  uranium. 
Comparison  of  the  spherical  shell  results  with  those  of  the  ECEL  test  sand¬ 
wiched  foils  is  limited  by  the  important  difference  in  geometrical  arrangement 
of  the  scattering/absorption  material.  For  elastic  scattering  by  an  idealized 
thin  shell  surrounding  a  point  detector,  the  neutrons  scattered  in  are  exactly 
compensated  by  those  scattered  out.  Thus,  although  the  shells  used  by  Davey 
and  Amundson  are  not  thin,  their  results  would  not  be  expected  to  show  as  large 
flux  depression  effects  from  elastic  scattering  as  would  the  ECEL  results  using 
plane  foil  sandwiches.  > 

Davey  and  Amundson's  shell  measurements  were  .also  made  in  a  fast  re-' 

t  • 

actor  somewhat  similar  to  the  ECEL.  These  shell  measurements  and  >th'e  ECEL 
test  were  both  conducted  in  a  core  region  already  containing  large  quant’ ties  of 
aluminum,  carbon,  stainless  steel,  uranium,  etc.  It  is  likely  that  measured 
flux  depressions  caused  by  materials  placed  adjacent  to  a  foil  packet  might  be*, 
more  pronounced  in  environments  not  containing  as  much  of  these  materials; 
that  is,  in  an  initially  unperturbed  spectrum, 

4.  Croaa  Section  Re-Evaluation 

In  he  course  of  the  data  analysis  from  the  six  tests  studied,  it  became 
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evident  that  effects  of  the  Fe  (n,p)Mn  foil  activity  measurements  were  sub¬ 
stantially  inconsistent  with  those  of  all  other  foils  used.  This  wa9  also  found  to 

58  58  32  32 

be  the  case  for  the  Ni  (n,  p)Co  and  S  (n,  p)P  reactions,  but  to  a  lesser  de- 

32  58 

gree.  It  was  further  found  that  the  S  and  Ni  results  were  mutually  incon- 

54 

sistent  independently  of  the  Fe  measurements.  Figures  12  and  14,  and  Ta¬ 
bles  XIV  and  XVIII,  which  show  the  WMGR  and  GMRW  results,  respectively, 
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are  cited  as  examples  of  the  kind  of  effect  observed,  The  large  peak  (in'  the 
V/MGR  result)  and  dip  (in  the  GMRW  result)  at  1.7  Mev  in  the  spectrum  obtained 
with  the  "old"  cross  section  evaluation  are  due  to  what  can  be  considered  an  in¬ 
appropriate  curve  drawn  through  available  differential  cross  section  data,  espe¬ 
cially  in  a  region  near  the  reaction  threshold,  where  the  curve  is  quite  steep. v 

A  comprehensive  investigation  of  these  inconsistencies,  involving  many 

computer  runs  not  reported  here,  resulted  in  a  substantial  re-evaluation  of  the 
54  54  _ 

Fe  (n,p)Mn  cross  section  data.  This  investigation  is  reported  in  refer¬ 
ence  36.  Figure  23  shows  the  "old"  evaluation  used  (solid  curve)  as  well  as  the 
new  one  obtained  from  this  comprehensive  cross  comparison  of  all  tests  in- 
volved  (dot -dash  curve).  The  difference  between  the  two  curves  does  not  ap¬ 
pear,  perhaps,  to  be  significant  near  the  threshold  because  of  the  low  magni¬ 
tude;  it  was  seen,  however,  to  be  much  more  significant  than  the  apparently 
larger  difference  between  about  6  and  13  Mev,  because  of  its  much  greater  ef¬ 
fect  on  the  sensitivity  curve  (product  of  cross  section  and  differential  flux). 

Figures  12  and  14  compare  the  solution  spectra  resulting  from  the  newer 
evaluation  (top  curve  in  each  figure)  with  the  original  results,  for  these  two  tests. 
In  both  cases,  the  newer  results  are  seen  to  be  much  more  reasonable, 
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The  newer  Fe  (n,p)Mn  evaluation  is  currently  used  in  SAND  II  (see 
Volume  III);  it  is  not  presented  as  being  correct,  but  only  less  incorrect  than 
the  "old"  one  (more  properly,  it  gives  more  nearly  consistent  results  with  dif¬ 
ferent  kinds  of  spectra).  This  reaction  certainly  still  requires  further  re- 
evaluation. 


*After  the  completion  of  this  report,  it  was  found  that  a  narrow  peak  had  been 
inadvertently  introduced  into  the  KAFL  evaluation  at  1.75  Mev  and  was  the 
cause  of  the  sharp  peak  and  dip  in  the  WMGR  and  GMRW  solution  spectra  re¬ 
spectively.  The  presence  of  this  peak  in  the  KAPL  cross  section,  however, 
would  not  be  expected  to  significantly  alter  the  conclusions  reached  about  the 
Fe54(n,p)  c  ross  section. 

>!< ’!< T h e  new  evaluation  is  a  modification  of  that  reported  in  reference  17,  based 
on  a  heavier  weighting  of  the  experimental  data  of  Salisburg  and  Chalmers 
(see  Figure  23).  It  should  be  noted  that  the  reference  17  evaluation  was  origi¬ 
nally  replaced  by  the  "old"  (KAPL)  evaluation  (solid  curve  in  Figure  23).  The 
KAPL  evaluation  was  based  on  a  Haus er -Feshbach  calculation  up  to  5.8  Mev 
normalized  to  fit  the  data  of  Carroll  and  Smith,  and  an  extrapolation  to  higher 
energies  through  14  Mev  data,  and  had  originally  been  thought  to  be  more  ac¬ 
curate  than  the  reference  17  evaluation. 
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SECTION  VI 

ERRORS  IN  FLUX  SPECTRAL  DEFINITION  " 

A.  Summary  Discussion 

With  substantial  errors  in  cross  section  and  measured  activities,  an  exact 

fkl  - 

solution  to  Eqs.  (38)  is  likely  to  yield  values  of  <p  (E)  that  contain  large  errors, 
because  the  experimental  input  errors  can  cause  ill-conditioning  in  the  system 
of  equations.  Gold  (Ref.  8),  DiCola  and  Rota  (Ref.  11),  and  McElroy  (Ref.  1) 
have  investigated  the  propagation  of  these  experimental  errors  and  their  effect 
on  the  accuracy  of  different  types  of  exact  solutions  for  neutron  spectra.  Grundl 
(Ref.  10),  DiCola  and  Rota  (Ref.  11),  and  Greer  and  Walker  (Ref.  4)  have  inves¬ 
tigated  methods  which  involve  approximate  solutions  for  the  neutron  energy  re¬ 
gion  above  a  few  hundredths  of  an  Mev,  while  McElroy,  Berg  and  Gigas  (Ref.  3) 
have  considered  the  same  problem,  but  for  the  entire  neutron  energy  region. 

Gold  (Ref.  9)  has  studied  the  problem  from  a  more  fundamental  point  of  view  by 
investigating  an  approximating,  but  general,  iterative  "unfolding"  method. 

These  latter  studies  have  shown  that  the  use  of  an  approximate  technique  for  ob¬ 
taining  an  appropriate  solution  to  the  system  of  equations  minimizes  the  effects 
of  the  propagation  of  experimental  errors.  The  SAND  II  iterative  procedure 
was  designed  to  accomplish  this  type  of  minimization. 

As  part  of  the  error  analysis  study,  the  SAND  II  code  was  used  to  analyze 
available  data  for  several  different  types  of  neutron  environment,  as  discussed 
in  Section  V.  The  results  of  these  tests  as  well,  as  the  analytical  tests  discussed 
in  Section  IV  indicate  that  the  code  is  capable  of  producing  a  solution  with  rea¬ 
sonably  correct  envelope  of  the  differential  structure.  The  accuracy  with  which 
this  can  be  done  depends  strongly  on  several  factors;  principal  among  these  are 

1)  adequacy  of  energy  range  coverage  by  the  foil  set  used, 

2)  accuracy  of  foil  activity  measurements  (including  self-absorption  and 
scattering  effects), 

3)  accuracy  of  the  differential  cross  section  evaluations  included  in  the 
cross  section  library,  and 

4)  aptness  of  initial  spectral  approximation. 
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In  ill-conditioned  cases  (in  the  sense  of  any  of  the  four  items  above),  the 
code  appears  to  introduce  somewhat  unreasonable  perturbations  in  the  flux  spec¬ 
trum  over  certain  energy  regions.  A  good  deal  of  cross  comparison  among  all 
the  tests  analyzed  showed  that  it  is  often  possible  to  distinguish  among  the  var¬ 
ious  possible  causes  of  the  ill-conditioning,  when  such  unreasonable  structure 
occurs.  For  example,  anomalous  structure  in  a  solution  spectrum  can  probably 
be  attributed  to  inadequate  foil  coverage  if  it  occurs  in  an  energy  region  not  well 
within  the  sensitivity  limits  of  more  than  one  foil  reaction  used;  this  can  be  de¬ 
termined  by  examining  the  list  of  energy  limits  within  which  90%  of  the  activity 
is  produced  for  each  foil  (included  in  the  printed  output  of  each  run).  Such  an 
examination  might  result  in  a  decision  to  irradiate  additional  foils,  to  help  fill 
in  the  energy  "gaps". 

An  ill-chosen  initial  approximation  can,  on  the  other  hand,  result  in  sensi¬ 
tivity  limits  for  certain  foils  which  are  very  much  displaced  in  energy  (for  the 
approximation  spectrum)  from  their  true  location;  this  can  cause  iterative  cor¬ 
rections  to  be  made  in  inappropriate  energy  regions  resulting,  in  turn,  in  ap¬ 
proach  to  an  incorrect,  although  arithmetically  reasonable,  solution.  This  ef¬ 
fect  is  especially  obvious  in  the  case  of  resonance  detectors;  if,  for  example, 
the  true  spectrum  includes  a  substantial  flux  in  the  region  of  a  major  resonance, 
and  the  initial  approximation  does  not,  the  calculated  sensitivity  limits  can  be 
shifted  upwards  as  much  as  several  decades  in  energy  to  a  region  of  very  small 
cross  section,  to  compensate  for  the  loss  of  calculated  activity  under  the  reso¬ 
nance.  Depending  on  the  particular  set  of  foils  used,  this  may  cause  foils  which 
in  fact  did  overlap  in  energy  regions  of  sensitivity  to  appear  to  be  affected  by 
separate  regions,  or  vice  versa;  this  can  sharply  change  the  correction  factor 
averaging  among  the  various  foils,  and  produce  artificial  spectral  structure. 

An  inaccurate  foil  activity  measurement  (outside  of  the  bounds  of  acceptable 
experimental  error  in  a  particular  case)  can  usually  be  identified  by  standard 
statistical  techniques,  as  a  result  of  examining  the  difference  between  calculated 
and  measured  activities  for  all  foils  used  for  the  solution  spectrum.  Experi¬ 
ence  to  date  suggests  that  errors  in  activity  measurements,  particularly  with 
threshold  detectors  (or,  in  general,  with  detectors  which  overlap  strongly  in 
energies  of  sensitivity),  usually  result  in  perturbations  much  less  severe  and 
broader  in  energy  than  those  caused  by,  say,  inaccuracies  in  differential  cross 
section. 
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Inaccurate  differential  cross  section  data  evaluation  for  a  particular  reac¬ 
tion  causes  unreasonable  structure  which  is  not  specific  to  any  one  test.  Eval¬ 
uation  of  a  large  body  of  data  can  identify  a  reaction  which  results  in  systemat¬ 
ically  perturbed  solution  spectra,  or  for  which  a  correlation  is  recognizable  be¬ 
tween  the  nature  of  the  perturbation  and,  perhaps,  the  kind  of  environment  (ther¬ 
mal,  fission,  fusion,  etc.  ).  The  appearance  of  such  unreasonable  perturbations 
in  the  spectrum  is,  then,  actually  an  important  advantage,  since  it  can  be  used 
to  adjust  currently  used  cross  section  evaluations  in  the  SAND  II  code  to  provide 
more  consistent  results  for  all  kinds  of  spectra. 

It  is  recognized  that  this  discussion  of  distinguishing  among  various  causes 
of  ill -conditioning  is  over-simplified  for  the  sake  of  illustration.  In  practice, 
combinations  of  difficulties  are  generally  in  effect,  and  require  more  careful 
and  studious  interpretation.  It  may  be  difficult  to  specify  the  error  in  a  solution 
spectrum  resulting  from  inadequate  foil  coverage  and/or  aptness  of  the  initial 
spectral  approximation.  In  general,  some  knowledge  about  the  anticipated  form 
of  the  environment  under  study  will  be  required  to  make  a  proper  assessment  of 
such  errors. 

Errors  associated  with  cross  sections  and  foil  activities  are  discussed,  in 
somewhat  more  detail,  in  the  following  two  sections. 

B.  Evaluated  Cross  Sections 

Studies  have  been  made  of  the  accuracy  of  foil  detector  cross  section  data 
by  Grundl  (Ref.  10),  Beauge  (Ref.  37),  Liskien  and  Paulsen  (Ref.  38),  Barrall 
and  McElroy  (Refs.  17,  39),  Parker  (Ref.  40),  Boldeman  (Ref.  12),  Zijp  (Refs. 

6,  7)  and  others.  A  review  of  the  results  of  the  more  recent  of  these  studies 
supports  the  conclusion  that  the  integrated  accuracy  of  present  evaluated  cross 
section  data  is  within  ±  10%  for  many  foil  detector  reactions  in  current  use,  par- 

i\< 

ticularly  resonance  and  fission  reactions.  It  is  much  more  difficult,  however, 
to  assign  limits  of  accuracy  to  the  differential  data;  accurate  differential  cross 
section  evaluations  are  essential  to  determining  neutron  flux  spectral  shape, 


*This  means  that  errors  in  the  cross  section  data  for  many  of  the  foil  detector 
reactions  used  do  not  result  in  more  than  ±  10%  errors  in  the  activation  inte¬ 
grals  (over  the  pertinent  energy  ranges)  for  a  neutron  spectrum  such  as  l/E 
or  fission.  (See  references  10,  17,  and  39.  ) 
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since  this  requires  independence  of  the  response  functions  for  the  various  foil 
reactions  used. 

As  indicated  in  Section  V-B,  erroneous  (or,  more  properly  for  these  pur¬ 
poses,  inconsistent)  cross  section  evaluations  can  lead  to  important  spectral 
structural  anomalies.  It  is  important,  therefore,  that  currently  used  cross 
section  evaluations  be  surveyed  to  determine  which  contain  errors  that  contrib¬ 
ute  to  erroneous  spectral  structure  in  iterative  solutions,  and  that  these  errors 
be  eliminated  by  further  studies  and/or  re-evaluation.  The  multiple  foil  iter¬ 
ative  technique  may  be  applied  as  mentioned  in  Section  V-B, "to  unfolding  the  re¬ 
lative  shape  and  absolute  magnitude  of  differential  cross  section  functions. 

Dierckx  (Ref.  41)  has  stated  that  an  error  of  ±  5%  to  ±  10%  in  the  measure¬ 
ment  of  cross  sections  is  sufficient  for  measuring  fluxes  to  ±  1%,  when  the 
measurement  is  done  by  comparing  the  unknown  spectrum  with  a  known  spec¬ 
trum  which  is  only  slightly  different.  Grundl  (Ref.  10)  has  used  this  compari¬ 
son  procedure  in  studying  fission  type  spectra,  and  suggests  that  spectral  anal¬ 
ysis  with  activation  detectors  always  involves  a  spectral  comparison.  The  iter¬ 
ative  perturbation  technique  is  essentially  a  comparison  procedure,  since  its 
initial  step  is  the  selection  of  an  input  spectral  approximation,  which  may  be 
based  on  theoretical  predictions.  Dierckx's  comment  may  therefore  be  taken 
to  apply  to  the  iterative  perturbation  technique  to  some  degree,  depending  on 
the  accuracy  of  the  initial  approximation. 

C.  Activity  Measurements 

Few  statements  on  the  accuracy  of  the  measurement  of  foil  detector  activ¬ 
ities  have  appeared  in  the  literature  (Ref.  24).  For  many  reactions  accuracies 
to  within  ±  5%  are  obtainable.  It  appears  that  for  other  reactions,  however,  an 
accuracy  to  within  ±  10%  is  more  realistic.  In  the  multiple  foil  irradiations 
discussed  in  Section  V,  such  accuracy  was  more  often  the  exception  than  the  rule. 

Measurements  of  activity  are  often  affected  by  such  factors  as  impurity  ma¬ 
terials  which  contribute  interfering  decay  products,  incomplete  knowledge  of  de¬ 
cay  schemes,  uncertainties  in  half-life  and  inconsistencies  in  half-life  values 
used  among  various  laboratories,  and  the  lack  of  readily  available  standards  for 
products  with  short  half-lives. 
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One  example  of  the  importance  of  establishing  better  information  on  foil 

measurement  reliability  is  given  in  the  results  of  the  W'MGK  test  (Fig.  1.1). 

Table  XIV  shows  the  measured  and  calculated  activity  results  for  this  test.  It 
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is  evident  that  the  magnitude  of  uncertainty  in  the  measured  Fe  (n,p)Mn  ac¬ 
tivity  is  of  the  same  order  as  the  magnitude  of  inconsistency  of  the  ratio  of  the 
measured  to  calculated  activities  for  that  foil  with  the  others.  Conclusions 
reached  about  the  particular  structure  introduced  in  the  solution  spectrum  at 
1.7  Mev  (regarding  the  inaccuracy  of  the  cross  section  evaluation  for  this  reac¬ 
tion)  would  become,  therefore,  questionable  to  the  degree  to  which  these  con¬ 
clusions  were  unsubstantiated  by  other  results.  Obviously,  a  more  careful  at¬ 
tempt  to  achieve  accuracy  in  activity  measurement  in  this  case  would  have  been 
extremely  valuable.  (It  should  be  pointed  out  that,  in  the  actual  execution  of 
this  test,  there  was  no  intent  to  use  the  measurements  as  they  have  subsequently 
been  used;  the  test  was,  rather,  initially  only  a  check  for  material  impurities  in 
the  foils*  ) 

It  is  possible  to  improve  the  reliability  of  activity  measurements  by 
careful  interlaboratory  comparison  and  cooperation.  The  round  robin  study  of 
three  reactions  conducted  by  the  ASTM  in  which  eight  laboratories  participated 
is  a  good  example  of  how  this  can  be  accomplished  (see  Table  XIII).  For  the 

CQ  C  Q  CA  C  A 

Ni°°(n,p)Co  0  and  Fe3  (n,p)Mn  reactions,  agreement  among  all  the  labora¬ 
tories  was  within  a  spread  of  4-1/2%  of  the  average;  for  the  Co^(n,y  )Co^ 
reaction,  the  spread  was  10%  of  the  average;  except  fer  one  laboratory,  how¬ 
ever,  this  reaction  also  had  a  4-1/2%  spread.  Improvements  in  foil  activity 
measurement  reliability  are  necessary  and,  as  this  round  robin  study  shows, 
are  possible. 


SECTION  VII 
CONCLUSIONS 

The  analytical  and  experimental  results  presented  in  this  report  have  pro¬ 
vided  an  indication  of  the  degree  to  which  the  SAND  II  multiple  foil  iterative 
method  may  be  useful  in  achieving  the  immediate  and  long-range  objectives  of 
this  work,  as  stated  in  Section  II-A.  The  experimental  tests  have  been  partic¬ 
ularly  valuable,  because  they  have  involved  different  laboratories  and  have 
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helped  to  establish  a  true  measure  of  the  current  state-of-the-art.  The  results 
of  these  tests,  when  viewed  as  a  uniform  hotly  of  information,  strongly  support 
the  anticipated  accuracy  of  the  multiple  foil  iterative  method.  The  method  can 
be  expected  to  give  integral  neutron  flux  results  over  the  energy  range  from 
IQ’ 10  to  18  Mev  that  are  accurate  to  within  10  to  30%  at  all  energies,  and  is 
rather  uniquely  suited  for  the  determination  of  total  absolute  integral  flux.  Dif¬ 
ferential  flux  spectra  are  also  obtained  which  reflect  true  structure  as  well  as 
fluctuations  resulting  from  errors  in  foil  reaction  cross  sections  and  activities. 
The  method  is  such  that  the  use  of  a  large  number  of  foils  with  overlapping  en¬ 
ergy  regions  of  sensitivity,  and  subsequent  examination  of  the  solution  differen¬ 
tial  spectral  structure,  can  be  used  to  help  distinguish  the  true  structure  from 
that  which  may  be  caused  by  cross  section,  activity,  and  foil  self  absorption/ 
scattering  errors. 

For  some  tests,  the  foil  activation  results  have  been  compared  with  time- 
of -flight  spectrometer  measurements  and  Monte  Carlo  calculations,  while  for 
others  the  comparison  has  been  with  reactor  physics  calculations,  using  diffu¬ 
sion  and  transport  theory.  These  comparisons  indicate  that  the  multiple  foil  in¬ 
tegral  flux  results  are  in  good  agreement  with  the  time -of -flight  measurements 
and  Monte  Carlo  calculations  at  most  energy  points,  and  that  some  of  the  dis¬ 
agreement  that  does  exist  can  be  attributed  to  positioning  and/or  neutron  ab¬ 
sorption  and  scattering  effects.  In  general,  the  agreement  with  reactor  physics 
calculations  is  not  as  good,  except  for  those  cases  in  which  the  "real"  environ¬ 
ment  can  be  represented  by  an  appropriate  mathematical  model.  It  is  concluded 
that  these  studies  have  yielded  quite  credible  integral  flux  results,  which  have 
increased  confidence  in  (1)  the  accuracy  of  current  cross  section  evaluations  and 
(2)  the  ability  of  the  iterative  method  to  be  used  to  determine  the  reliability  of 
theoretical  predictions. 

The  addition  of  gold  as  the  third  cover  material  in  the  SAND  II  code  has 
permitted  a  partial  study  of  the  effectiveness  of  using  self-sandwiched  resonance 
detectors,  by  allowing  calculations  involving  self- sandwiched  gold.  Results  in¬ 
dicate  that  this  stratagem  is  quite  practicable;  the  effect  is  to  alter  the  response 
of  the  resonance  detectors  so  that  their  behavior  is  more  nearly  like  that  of  in¬ 
tegrating  threshold  detectors.  The  self-sandwiching  technique  is  limited,  how¬ 
ever,  by  the  approximations  made  in  calculating  attenuation  by  the  sandwich 
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material;  the  degree  to  which  these  approximations  are  acceptable  depends 
strongly  on  the  geometry  of  a  particular  case.  Careful  study  must  be  made  of 
the  characteristics  of  the  materials  used  as  resonance  detectors  in  the  partic¬ 
ular  environment  in  which  they  are  used. 

The  authors  have  concluded  from  results  of  analyzing  spectra  that,  in  addi¬ 
tion  to  providing  experimental  measurement  of  neutron  flux  spectra,  the  iterative 
method  employed  affords  a  procedure  for  experimental  evaluation  of  foil  reaction 

differential  cross  section  data.  Studies  are  reported  in  which  the  authors  used 
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the  code  to  re-evaluate  the  differential  cross  section  data  for  the  Fe  (n,  p)Mn 
reaction.  (Because  of  its  long  product  half-life,  this  reaction  is  often  used  for 
long-term  irradiations  to  measure  fluence  [time-integrated  neutron  exposure] 
and  is,  therefore,  of  considerable  theoretical  as  well  as  practical  interest.  )  It 
is  not  suggested  that  the  re-evaluation  is  correct  in  either  absolute  magnitude  or 
differential  form;  it  is,  after  all,  without  detailed  structure,  and  somewhat  arbi¬ 
trarily  drawn  in  several  energy  regions.  This  evaluation  is  offered  as  the  best 
currently  available,  only  in  the  sense  that  it  seems  to  produce  the  most  nearly 
self-consistent  results  in  the  tests  to  which  it  has  been  applied. 

Because  of  the  presence  of  structural  materials  in  a  neutron  environment, 
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flux  depressions  can  exist  which  affect  such  reactions  as  Fe  °(n,  p)Mn  , 

54  54  58  58 

Fe  (n,  p)Mn  ,  and  Ni  (n,  p)Co  ;  results  reported  here  indicate  that  current 
differential  cross  section  evaluations  for  these  and  other  reactions  are  not  suf¬ 
ficiently  detailed  to  allow  the  SAND  II  code  to  "unfold"  these  depressions  prop¬ 
erly,  and  require  more  detailed  study. 
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